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ABSTRACT 


\ 


This  ri'port  prosonts  _l^survoy  of  th(“  various  phases  of  calculations  which  led  to  the  quan¬ 
titative  prediction  of  the  important  undei  water  explosion  parameters  for  Operation  Wigwam.  i 

Vit  I'.iiiumnn  of-^ttrrlrfirr  Water.  The  analysis  of  the  explosion  phenomena  requires  a 
k.iowledge  of  the  thermodynamic  properties  of  water  over  an  extremely  wide  pressure  range, 
i.e.,  from  infinity  down  to  the  low  pressures  of  an  acoustic  wave.  No  equation  ol  state  is 
Known  which  satisfactorily  covers  this  range;  therefore  five  separate  pressure  ranges  were 
ronsidercxlk 

Region^.  At  eslramely  nlgh-prasturcs  and  tempeiatwecii  the  moleciilca. of-waUr  ar«.CMn- 
pletely  diss  ciated  and  ionized.  The  gas  is  ideal  and  n^atomic  if  tKe  small  effects  orradia- 
tion  pressut »  and  electrostatic  forces  are  excluded.^'  \  / 

Region  1.  For  somewhat  lower  pressures  and^mperatures  the  ^edium  is  only  parj^ally 
dissociated  ind  ionized.  LaboiVous  equilibrium  emulations  were  made  to  determine  the! 
thormodyn/iiiic  data  in  this  region.  The  p-v-T  r^Ution  necessary  for  tills  purpose  was  oji- 
tained  fryfc  the  detonation  theory  of  high  cxplos^es,  In  particular  from  the  explosive  hyirazlne 
nit  rat vyAich  forms  H-0  as  its  wincipal  reactmn  product.  \  / 

Jj^.on  ill.  At  still  lower  pressures  anS  btmperatures  the  water  molecihe  remalnsylntact. 
Fy^hii  range,  calculations  usinir.the  ThonjM-Fermi-Dlrac  theory  were  mat^  / 
y'  i.utioiis  IV  .Hid  V.  For  pressurhsfiyo  723,000  psl  down  to  acoustic  values/^ttmt  ex- 
/ ^lerimnitai  measurements  notably  by'Kl^mtn*  and  Ctrnevale  find  Lltovltz**  were  used. 

^  Shoe':  Ware  I’ltciiomcmi.  Tliiisw  i  iliiiil  iliiinw  aeif  iliiii  srjiaiaiail  lulu  ||i'v‘‘’~at 
'**ror  extremely  high  pressures  the  solution  of  the  point  blast  problem  of  Taylor^^^  applicable. 
For  lower  pressures  the  three  partial  differential  equations  of  the  spherical  fluid  motion  were 
Integrated.  The  method  was  not  tractable  below  a  shock  pressure  of  about  450,000  psl  (corre¬ 
sponding  to  a  shock  radius  of  81  ft  In  Operation  Wlgwamh 

low  prcssii^sJ)SJuaa«*-of  theSllay-Wallhras’  sHock-wJe  theory.  At  very  low  pressures 
asyiTtptCric  relations;  similar  to  those  first  declycf^JigrXlrkwood  aitd  Beth*,*  w«re  used.  ~ 


limilar 

BttbtiU.£lu:iioiK*:ua.^The  energy  dissipation  I 


1  at  the  frtmt  .if  the  intense  shock  wave  from  a  point  explosion  produces  the  heat 
wiiicii  vaporizes  the  water  and  foi  iiis  a  Steam-filled  cavity.  This  bubble  pulsates  in  a  manner 
similar  to  that  observed  for  bubbles  produced  by  high  explosives.  The  analysis  yielded  the 
maximum  bubble  radius  and  the  period  of  the  first  pulsation  as  well  as  the  total  mass  of  water 


evaporated  up  to  the  moment  of  the  first  bubble  maxlmunr 
clwling  the  rapid  up^Ktrd  migratioitrcan  be.xaii 
This  establishes  ax  upper  limit  for  the  periods  and  ter 
actual  behavior^ stvam  lHrl,bles  has  been  studieir-witll^ 
energy  sourc^^  The  results  of  these  >ests  wej*  used  to 
condensatlnn.'Vihich  occurred  in  Wigwam.  It  tinned  out  that 
been  condensed  before  the  bubble  reached  tKe  surface  and  thai 
some  ry.s^hlanee  to  the  "breakthnggjlill 


lOdf 
I  obti 


>  later  buUblu  {JllinTomenaJnr* 
ior  high-exploslye'gas  bubbles, 
ligratlon  ol  a  steam  bubble.  The 
tests  using  electric  sparks  as 
Informatloti  on  the  amount  of  . 
ost  all  the  vapor  must  have  ! 
surface  phenomena  which  lad 


IS^  by-U>e-viotent  -upwellftig 


j 

1 


o{4h£-aiat«'  which  previously  surrounded  the  bubble  and  which  acquired  the  latter’s  upward 
moinentum. 

In  sumniat  y,  it  was  found  that,  In  the  region  where  pressures  are  less  than  3000  psl,  the 
calculated  pressure-d'stance  curve  is  similar  to  one  from  TNT  ha^‘  x  about  69  per  cent  as 
much  energy.  In  this  stmc  region  me  calculated  shock-wave  energy  flux-difa.-uicc  curve  is 
similar  to  one  from  TNT  having  about  82  per  cent  as  much  energy.  The  maximum  bubble 
radius  was  calculated  to  be  3 16  ft  with  a  first  bubble  period  of  2.88  sec.  This  period  corre¬ 
sponds  to  that  from  TNT  .'laving  8i  per  cent  as  much  energy.  The  amount  of  water  evaporated 


was  calculated  to  occupy  the  same  volume  as  a  30-kt  TNT  s 
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PREFACE 


Project  1.1  of  Operation  Wigwam  was  one  of  four  projects  (i.l,  1.2,  1.4,  and  1.5)  lor 
which  the  Nasal  Ordnance  Laboratory  was  responsible.  Its  objective  was  to  determine  tlie 
principal  underwater  explosion  phenomena  to  be  expected  from  the  explosion,  at  a  depth  of 
2000  ft  in  deep  water,  of  an  atomic  device  having  a  nominal  yield  of  30  metric  kllotons.  The 
successful  achievement  of  this  objective  enabled  at  least  two  extremely  Important  practical 
results  to  become  available:  (1)  the  determination  of  proper  locations  for  the  targets  and 
instrumentation  during  the  Operation  and  (2)  the  development  of  methods  for  predicting  under¬ 
water  explosion  phenomena  from  other  yields  and  firing  geometries. 

This  summary  report  not  only  gives  the  predictions  which  were  used  in  helping  to  deter¬ 
mine  the  experimental  configuration  but  makes  comparisons  of  these  predictions  with  the  actual 
measurements  obtained.  The  good  agreement  between  theory  end  e.\perlment  lr.dicates  that  the 
methods  used  describe  the  important  phenomena  with  satisfactory  accuracy. 

In  this  report  the  important  equations  which  have  been  used  in  Project  i.i  are  summariaed 
and  explained.  Only  simple  derivations  are  given.  For  a  complete  aiulysis  the  following  re- 
port.s,  which  describe  the  subject  matter  more  thoroughly,  should  be  consulted: 

NAVOKO  Report  4181:  An  Equation  of  State  for  Water,  by  Hans  G.  Snay  and  John  F.  fiutler 
(in  preparation). 

NAVORD  Report  3847:  An  Equation  of  State  lor  Water  at  Extreme  Pressures,  by  J,  H. 
Rosenbaum. 

NAVORD  Report  4182:  A  Theory  of  the  Shock  Wave  Produced  by  a  Point  Explosion,  by 
Hans  G.  Snay  (in  preparation). 

NAVORD  Report  4183:  An  Analysis  of  Solutions  of  the  Point  Blast  Problem,  by  Andrg  N. 
Glcyzal  (in  preparation). 

NAVORD  Report  4184:  Numerical  Analysis  of  the  Underwater  Point  Blast  Problem,  by 
John  F.  Butler  (in  preparation). 

NAVORD  Report  4185:  Underwater  Explosion  Phenomena  U:  The  Parameters  of  Migrating 
Gas  Bubbles,  by  Hans  G.  Snay  (in  preparation). 

The  reader  who  is  not  interested  in  mathematical  details  is  invited  to  read  the  introduc¬ 
tory  and  summary  paragraphs  of  each  chapter  of  this  'eport,  as  well  as  Secs.  3.3  through  3.5. 
riiis,  together  with  a  study  of  7'abie  3.1.  the  figures,  and  the  glossary  (Appendix  A),  will  pro¬ 
vide  a  fair  idea  of  the  methods  used  and  the  results  obtained  in  this  project. 
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CHAPTER  1 

EQUATION  OF  STATE  FOR  WATER 


1.1  INTRODUCTION 

An  annlyHlB  nf  thn  pxploslon  phenomenn  following  the  firing  of  an  atomic  device  under 
wntoi'  requircK  the  knowledge  of  the  thermodynamic  properties  of  water  over  a  range  extending 
(roni  the  exlremcly  high  presBuren  and  temperatures  occurring  Immediately  after  the  explosion 
d'i'vn  to  the  eomititonB  exinilng  in  an  acousllr.  wave.  No  satisfactory  equation  of  slate  is 
kno  >/M  whieh  Rovers  the  whole  range  of  interest. 

In  the  Initial  phase,  subsequently  called  Region  I,  the  temperature  is  so  high  that  the  atoms 
at  e  (  oniplctely  stripped  of  all  their  electrons. 

At  somewhat  lower  temperatures  (Region  II),  the  medium  consists  of  a  mixture  of  the 
foll'iwlng: 

1.  Diatomic  molecules  and  radicals  formed  from  hydrogen  and  oxygen  (OH,  Hj,  O],  etc.). 

2.  Monatomic  hydrogen  and  oxygen.  , 

3.  Hydrogen  and  oxygen  Ions  (O*,  H",  etc.). 

■I,  ?'ree  oiccirons. 

i),  lonixed  molecules  and  radicals  (Oj,  OH",  etc.). 

The  calculation  of  the  thermodynamic  properties  of  such  a  mixture  requires  a  knowledge 
of  the  concentration  of  the  various  constituents  of  the  mixture.  The  Halford-Klstlakowsky- 
VVilRon  (HKW)  equation  of  state  is  used  In  this  region. 

At  still  lower  pressures  and  temperatures  the  water  molecule  remains  Intact  (Region  III), 
hut  the  pressures  arc  still  far  above  the  range  where  direct  experimental  measurements  arc 
possible.  Ill  this  region  the  Thomas* Ferml'DIrac  (TFD)  theory  may  be  used.  However,  this 
theory  yields  acceptable  results  only  at  the  high>pressure  end  of  this  region.  Therefore,  to 
obtain  data  for  Region  III,  interpolations  must  be  made  between  these  calculated  high  pressures 
and  Region  IV,  for  which  measurements  by  Bridgman  are  available  (up  to  72B,000  psl).  Since 
the  Ranldnc-Hugonlot  curve  has  little  curvature  in  a  In  p  -  In  v  plot,  the  shock-front  data  can  be 
readily  interpolated  graphically. 

Region  IV  has  been  treated  In  several  publications.  Thermodynamic  data  behind  the  front 
are  obtained  by  using  a  modified  form  of  Ihe  Isentroplc  Talt  equation  and  adjusting  the  constant! 
In  this  equation  in  such  fashion  that  the  isentroplcs  fit  the  data  at  the  shock  front  (Rankine- 
Hugunlut  curve)  and  at  the  saturation  line. 

In  the  region  of  relatively  low  pressures  (Region  V),  where  the  shock  wave  behaves  almost 
like  .an  acoustic  wave,  the  thermodynamic  properties  of  water  may  be  Inferred  from  experi¬ 
ments  on  the  velocity  of  sound  In  water  ns  a  function  of  pressure. 
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1.2  REQUIRED  DATA 


For  any  bhork-wave  calculation  the  Ranklne-Hugonlot  parameters  must  bo  known  for  the 
medium  In  (4uestlon.  The  Ranklne-Hugonlot  conditions  are: 

the  Ranklne-Hugonlot  adiabatic 

E,-E:  =  5iJ--?£«(v,-v,),  (1,1) 

the  propagation  velocity  U 

pno  nartlri®  vploclly  u« 

(u, -u,)*^p,(v,-v,).  (1.3) 

The  subscript  1  designates  the  state  directly  behind  the  shock  front,  and  the  subscript  0  refers 
to  the  state  ahead  of  the  front.  E  is  the  internal  energy  per  unit  mass,  p  is  the  excess  pressure 
above  the  static  pressure  Pg,  and  v,  the  specific  volume,  is  the  reciprocal  of  the  density  p. 

The  Ranklne-Hugonlot  adiabatic  gives  the  p-v  relation  for  the  thermodynamic  change  of 
state  at  the  shock  front.  To  evaluate  this,  one  must  know  the  interrelation  between  internal 
energy,  pressure,  and  volume  A  formal  simplification  can  be  made  by  Introduction  of  the 
"reduced  Internal  energy” 

J  ’  J(p,v)  = (1.4) 

This  term  will  be  frequently  used  in  our  calculations.  It  is  a  dimensionless  magnitude  related 
to  the  heat  capacity.  For  an  ideal  gas  at  high  temperature  the  following  simple  equation  holds; 

where  the  superscript  0  indicates  the  Ideal-gas  state,  Cy  Is  the  beat  capacity  at  constant  volume, 
R  is  the  fi<Aa  constant,  and  r’  is  the  ratio  of  the  heat  capacities  at  constant  pressure  and  at 
constant  volume. 

The  Rankine-Hugar.iot  adiabatic  is  generally  giver  by 


Alter  rearrangement,  we  obtain  with  the  use  of  Eq.  1.5 

Pi  ' 

7TI ''I-''" 


This  simple  relation  holds  for  Ideal  gases  only.  For  real  gases  J  must  be  evaluated  from  the 
internal  energy. 


E  -  Eg  =  2  ntcl ,  dT  +  f  |t(|^)  -  (P  +  P,)  dv. 
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where  m  Is  the  number  ol  moles  of  the  Uh  constituent  In  the  mediuni,  Cv,  Is  Ite  Ideal'gas  best 
capacity,  and  T  is  the  absolute  temperature.  The  subscript  v  on  the  partial  differential 
quotient  in  the  second  integral  means  “at  constant  volume."  A  corresponding  notation  referring, 
for  instance,  to  constant  temperature,  entropy,  etc..  Is  used  in  this  report.  In  order  to  evaluate 
the  first  integral,  the  composition  of  the  medium  must  be  known.  This  information  will  be  ob¬ 
tained  from  the  equilibrium  calculations.  The  second  Integral  accounts  for  tlie  imperfect-gas 
behavior.  For  its  evaluation,  we  need  an  equation  of  state  for  the  imperfect  medium,  The  lower 
limit  of  this  integral  refers  to  the  specific  volume  at  the  ideal-gas  otate. 

In  considering  a  nonldeal  medium,  attention  must  be  given  to  the  use  of  the  symbol  y.  The 
customary  definition  of  this  symbol  is  y  =  Cp/Cy.  However,  in  the  hydrodynamic  literature  y 
stands  for  the  logarithmic  slope  of  an  Isentroplc.  in  general,  the  specific  heat  ratio  is  equal  to 
this  logarithmic  slope  only  for  the  case  of  an  ideal  gas.  In  this  report  the  symbol  y  will  mean 
the  logarithmic  slope  of  an  isentroplc,  l.e.. 


(1.9) 


where  c  sound  velocity  and  S  -  entropy. 

The  genera!  expression  for  this  isentroplc  exponent*  is 


y  -  - 


/3_{np\  ^  /a  In  p\^  /  /8E\ 

\3\nvl^  T  Vainv/T/  Wy’ 


(1.10) 


The  isentroplc  exponent  y  is  related  to  the  reduced  Internal  energy  J  by  the  following  simple 
equation: 


y 


1  P  '  Pq 
j  p 


+  1  + 


(1.11) 


For  c.inslant  J  this  equation  is  equivalent  to  the  ideal-gas  relation  1.5.  However,  in  the  general 
case  of  an  imperfect  gas,  Eq.  1.11  is  of  little  help.  It  is  simpler  to  calculate  y  and  J  inde¬ 
pendently  using  Eqs.  1.10  and  1.6. 

The  behavior  of  the  Ranklne-Hugonlot  adiabatic  can  be  expressed  in  the  follcwing  concise 
form: 


d  In  P|  _  d  In  Pi 
^  OH  d  In  V|  d  In  pj 


(1.12) 


This  definition  is  analogous  to  that  of  the  isentroplc  exponent  y  except  that  the  differential 
quotient  Is  taken  along  the  Ranklne-Hugonlot  adiabatic  Instead  of  the  isentroplc.  This  magnitude 
will  be  widely  used  In  the  hydrudynamlu  equations  of  this  report. 

The  equation  for  the  entropy  increment,  which  is  also  needed  in  the  hydrodynamic  calcu¬ 
lations,  Is 


The  dissipated  enthalpy  Increment  h  is  dsflned  by 
h  =  J  (Vo  +  V.)  -  J‘’‘  v{p,S  =  S,)  dp. 


(1.13) 


(1.14) 


‘The  word  "exponent"  Is  used  loosely  here.  For  constant  y,  Eq.  1.9  can  be  Integrated,  and 
it  yields  the  familial'  relation  pvi'-  const.  This,  however,  does  not  hold  for  a  variable  y. 
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rhi--  ma);iututie  ulves  the  lncreas,e  of  enthalpy  after  the  passage  of  the  shock  front  when  the 
medium  has  lseiitroptcal!\  returned  to  the  initial  pressure  Pg.  Another  expression  for  h  is 


h=  T(S,p-=0)dS 
=  J7,To(e‘^i'^-  1), 

where  Tj  is  the  temperature  of  the  water  before  the  explosion  and 

It,  cp  dT  _  .  1  r 

P  "  ■  T  -  To  ■  in  T/To  Jt.  T  ’ 

for  water  differ  only  very  slightly  from  the  value  of  the  actual  heat  capacity  Cp, 


(1.15) 


(1.16) 


1.3  REGION  OF  EXTREMELY  HIGH  TEMPERATURES  AND  PRESSURES  (REGION  1) 


The  equation  of  the  Rankiiie-Hugonlot  adiabatic  for  an  ideal  gas,  Eq.  1.7,  shows  that  the 
pri-"!siirp  becomes  Infinite  for  a  finite  value  of  v,.  Although  this  relation  holds  for  ideal  gases 
only,  it  is  generally  true  for  shock  waves  in  any  medium  that  V|  remains  finite  when  pressure 
and  temperature  approach  Infinitely  high  values.  At  such  conditions,  namely,  moderate  den¬ 
sities  but  infinitely  high  temperatures,  the  medium  is  completely  dissociated  and  ionized,  i.e., 
the  atoms  are  strippe.1  of  all  their  electrons.  The  medium  consists  only  of  such  small  particles 
as  electrons,  protons,  and  nuclei.  Each  of  these  has  but  three  degrees  of  freedom;  therefore 
I  =  *j.  The  small  particle  size  precludes  gas-lmperfectlon  effects,  and  such  a  plasma  would 
behave  like  an  ideal  gas  save  for  the  effects  of  radiation  pressure  and  electrostatic  lorces.  For 
our  purposes  it  seemed  permissible  to  neglect  these  effects  in  this  region  and  to  assume  that 
the  medium  behaves  like  an  ideal  monatomic  gas.  From  Eq.  1.7  wc  obtain 


Hill  v,  =  V( 
I'l-*-’ 


1  _  n 

r'  +  l  ‘  4  ■ 


The  shock-front  data  calculated  by  the  methods  described  in  the  following  paragraphs  were 
extended  to  infinite  pressures  by  considering  the  trend  of  as  a  function  of  v,.  The  magnitude 

1  _  d  In  V| 

■  “  d  In  ?, 

vanishes  at  v,  -  vg'4  according  to  the  above-discussed  behavior  ol  V|.  The  simple  technique 
used  is  described  in  ”AVORD  Report  4*81. 


1 .4  EQUATION  OF  STATE  IN  REGION  II 

To  evaluate  the  thermodynamic  functions  discussed  in  Sec.  1.2,  we  need  a  p-v-T  relation 
from  which  the  necessary  differential  quotients  may  be  determined.  No  completely  satisfactory 
equation  of  state  is  known  (or  media  in  Region  II,  where  dissociation  and  ionization  occur. 
Closest  to  this  range  are  the  thermodynamic  states  of  detonating  high  explosives  whose  reaction 
products  attain  pressures  up  to  4.3  million  psi,  temperatures  up  to  30fl0°K,  and  densities  above 
2  g  cc.  Several  attempts  have  been  made  to  describe  such  thermodynamic  conditions  by  means 
cf  an  equation  of  state.  Ail  these  equations  have  objectionable  features  from  a  theoretical  point 
of  view.  However,  the  theory  of  the  detonation  process  permits  a  determination  of  the  arbitrary 
constants  in  such  equations  from  experimentally  measured  detonation  rates.  These  rates  in¬ 
crease  with  Increasing  loading  density  of  the  explosive,  an  effect  which  is  solely  due  to  the 
imperfect-gas  behavior  of  the  reaction  products.  Therefore  this  is  a  sensitive  method  of  de¬ 
termining  imperfection  teims  of  high-pressure  equations  of  state. 
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For  this  project  we  have  chosen  the  HKW  equation  ot  state  because  this  equation  has  been 
tcstc'i  aitalnst  <!elonatloii  data  o(  a  (treat  number  of  explosives  with  good  success.*  Ip  particular, 
thi.s  equation  predicts  the  detonation  rate  of  the  explosive  hydrazine  nitrate  with  satisfactory 
accuracy.  The  dominant  reaction  product  of  this  explosive  Is  HjO.  Hence  one  should  expect 
that  this  equation  is  applicable  to  onr  problem  within  the  limits  of  Its  validity. 

The  »IKW  equation  of  state  has  the  form: 


(p  +  Po)v  2^1  (1  1  .-ce***) 


k  S  n.kj  ng  =  £  nj  «  =  0.25  =  0.3, 

I  : 

where  n,  is  the  number  of  moles  of  the  ith  component  in  Mg  weight  units  of  gas  and  kj  is  Its 
covulunic  factor.  A  few  covolume  factors  have  been  determined  from  the  measured  deto.iation 
rates.  For  our  problem  it  was  necessary  to  estimate  the  additional  covolume  factors  which 
could  no*  he  determined  by  this  method.  A  complete  account  of  this  is  given  In  NAVORD  Report 
4181,  where  a  list  of  the  numerical  values  is  found.  A  few  of  the  more  Important  components 
arc  quoted  here: 


Component 

kl 

HjO 

285 

O2 

300 

Hj 

60 

ot 

255 

OK 

200 

Component 

k| 

0 

100 

H 

20 

0+ 

80 

65 

H+ 

0 

The  limits  of  validity  of  the  HKW  equation  are  given  by  the  Impeifcction  term  x.  Good 
rcsult.s  were  obtained  In  the  study  of  detonation  phenomena  for  valves  of  x  between  1.7  and  4.2. 
This  range  is  determined  by  the  experimental  data  available.  It  Is  not  unreasonable  to  assume 
that  acceptable  accuracy  may  be  obtained  for  much  lower  values  of  x,  since  Eq.  1.17  reduces 
to  the  perfect-gas  law  foi  small  x. 

If,  within  these  limits  of  x,  the  HKW  equation  Is  a,,plled  to  water,  higher  pressures  and 
tenipcrnturo.s  result  than  those  usually  obtained  for  explosion  reaction  products  witi.  the  same 
X.  This  is  due  to  the  low  molecular  weight  of  water  and  the  relatively  high  number  of  particles 
in  the  dissociated  state. 

The  covolumc  (actors  noted  above  hold  for  Ihe  comparatively  low  temperatures  occurring 
i,i  the  detonation  of  high  explosives.  One  may  saivly  assume  that  the  constituent  molecules  and 
atoms  are  In  the  jtround  electronic  state  under  these  conditions.  This  raises  the  question  of 
whether  these  covolume  factors  arc  applicable  to  the  much  higher  temperatures  occiuring  in 
Region  II,  where  appreciable  excitation  and  Ionization  are  encountered. 

It  is  well  known  that  the  partition  function  for  an  Ideal  gas  diverges  at  high  temperatures. 
However,  it  Is  possible  to  show,’’^  that  convergent  series  for  the  partition  function  are  obtained 
for  imperfect  gases.  If  we  consider  a  HKW  gas  consisting  ol  only  one  component  and  regard  the 
various  excited  states  as  separate  species,  the  partition  function  is  given  by 


0  ^ 


P 

1  -f 


tii  exp 


(1.18) 


where  k'  =  Boltzmann  coiistant 
h  Planck’s  constant 
ni  -  particle  mass 

gi  statistical  weight  ol  the  llh  excited  state  (l.e.,  of  the  Uh  species  of  the  mixture) 
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{|  •  excitation  encrjty  ot  the  1/A  excited  state 

k*  covolume  lactur  for  the  ith  excited  state 

n*  =  number  of  moles  ot  atoms  which  are  In  the  t/A  excited  state 

n*  =  ::  ni  =  total  number  uf  moles  of  particles  in  M«  grams  of  the  mixture 

k*  ■•=  nt  kj  =  total  covulume  factor  for  the  mixture. 

(The  asterisks  are  used  here  to  emphasise  the  fact  that  we  are  now  considering  a  gas  consist¬ 
ing  cf  only  one  component,  for  Instance  O^,  and  are  regarding  the  various  excited  states  of  this 
component  as  separate  species).  For  x  -  0,  Eq.  I.IB  reduces  to  the  familiar  partition  function 
of  an  Ideal  gas.  In  this  •■ase,  the  sum  In  Eq.  1.18  Is  replaced  by  ^  gi  exp  (<i/k'T),  which  Is 
divergent  since,  although  the  Boltzmann  factors  exp  (<i/k'T)  become  smaller  with  Increasing  1, 
the  statistic  al  weights  gi  become  arbitrarily  large.  For  the  HKW  gas,  however,  the  last  factor 
in  the  sum  in  Eq.  l.i8  ensures  rapid  convergence. 

To  determine  how  many  terms  are  required  in  the  evaiitatlon  ot  the  partition  function. 

Eq.  1.18,  we  assume  that  the  covolume  factor  for  a  given  excited  state  of  an  atom  or  Ion  Is 
propo  "ir  --1  to  ths  effective  atomic  or  Ionic  volume  cl  the  atom  or  Ion  when  in  this  excited 
si;  .  'h'n  yjssible  (NAVORD  Report  4181)  to  derive  the  iollowlng  expression: 


k* = y  nui = kjnj  ^ - 1 , 


where  Vg  =  effective  acomic  or  ionic  volume  of  the  component  when  In  Its  ground  state 

Vi  =  effective  atomic  or  Ionic  volume  of  the  component  when  In  its  i/A  excited  state 
‘if  -  covolume  factor  for  the  component  when  In  Its  ground  state. 

An  evaluation  of  this  equation  yielded  the  interesting  result  that,  lor  all  conditions  where  the 
HKW  equation  ot  state  is  applicable,  only  those  excited  states  which  have  effective  atomic 
volumes  equal  to  or  near  that  of  the  ground  state  contribute  significantly  to  the  partition  func¬ 
tion.  Higher  terms  were  found  to  be  entirely  negligible,  even  for  a  value  of  x  as  low  as  C.2. 

The  Ideal-gan  nronertles  'cr  the  high  temperatures  needed  In  the  calculation  of  the  Internal 
enerev,  the  adiabatic  exponent,  and  the  entropy  were  calculated  from  statistical  thermodynamics 
on  this  basis.  A  co.mpiei"  des^rlptiun  is  given  in  NAVORO  Report  4181,  where  extensive  tables 
cf  the  pertinent  thermodynamic  data  can  also  be  found. 

With  the  use  of  the  HKW  equation  ol  state,  the  thermodynamic  functions  of  Interest  to  us 
take  the  foiiuwing  form: 


1  11 
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Thp  change  of  nj;  and  U  with  temperature  T  or  with  the  specific  volume  v — as  well  as  the  con¬ 
centrations  nt  of  the  various  constituents — are  deteruiined  by  equilibrium  calculations. 


1.5  EQUILIBRIUM  CALCULATIONS 

The  equilibrium  constant  for  an  HKW  gas  is 


K, 

1 

where 

2.23  X  10*  /. 
^  "  v'i(DT)*»  V 

or 

^  _  2.589  X  10*^ 

■■  Unsold  correction, 


where  Kpj  =  Ideal-gas  equilibrium  constant  lor  the  J/A  reaction 
D  =  dielectric  constant 
Zj  =  valency  of  the  1/A  component 

Zj  “  valency  of  the  j/A  reaction  which  must  be  of  the  type  x  =  x^j  +  Zje  (e  denotes  the 
free  electron). 

The  Vi’s  are  the  coefficients  of  the  1/A  component  in  the  J/A  reaction  equation;  those  on  the 
right-hand  side  are  counted  positive,  and  those  on  the  left-hand  side,  negative.  The  sums  con¬ 
taining  vj  run  over  the  terms  of  the  reaction  equation.  To  Illustrate  this,  assume  the  following 
reaction: 


0  =  0**  +  2e. 


Here 

vo*i  =  +1 

Vp  =  +2 

E  vj  =  +2 
I 

and 


nnr‘=!i2rjLe. 


where  no+*  Is  the  number  of  moles  of  O**  ions,  n*  is  the  number  of  moles  of  electrons,  and 
no  Is  the  number  of  moles  of  oxygen  atoms  in  Mg  grams  of  the  mixture.  Zi  for  0**  is  2,  and 
lor  e,  -1.  Hence  2  vjZf  =  6.  Zj  of  this  reaction  is  2. 

The  total  number  of  moles  of  all  hydrogen  nuclei  in  ora  mole  of  dissociated  and  ionised 
water  must  be  equal  to  2,  and  that  of  all  oxygen  nuclei  must  be  equal  to  unity.  This  can  be 
Witten  ae  follows; 


y,  +  ^  r,ini  =  2 

yj+Erjini  -  1, 

'  17 
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'vh.Tc  r,i  IS  th'*  iiuiiil'ei  I't  hy(it<vcii  nuclei  in  the  ith  compound,  Cjj  Is  ihe  number  of  oxygen 
nuclei  .•  ,  for  OH.  r,  1  an-l  1),  and  where. 

y,  I  ar.-.'.ier  of  moles  of  hydrogen  atoms  In  Mj  grams  of  the  mixture 

>j  -luiiilier  of  .iiolf-s  of  un-lomzed  oxygen  atoms  in  Mj  grams  of  the 
mixture 

These  two  magnitudes  are  >  efercr.tc  values.  They  are  by  nature  the  same  as  certain  ni’s  for 
which  ri  IS  unity,  and  they  serve  as  Independent  variables.  A  third  reference  value  of  this  kind 
is 


yj  ’  niimbei  of  moles  of  electrons  In  grams  of  the  mixture  =  ^  njZi  (1.26) 
Then  the  conrcntration  of  the  1/A  compound  Is 

'>1  Kiy['>  y^'*  y?-.  (1.27) 

There  are  as  many  siirh  equations  as  there  are  reaction  equations.  For  our  calculations,  we 


hax 

•e  considered  the  followin'’ 

reactions 

} 

Ilcaction 

En 

E  ViKi 

1 

E  viZj 

S'" 

>■11 

fjl 

Zi 

4 

1 

0  -  O*  ’  e 

1 

-?0 

2 

1 

0 

1 

1 

2 

O  -  O*’  .  2e 

2 

-35 

6 

2.5876 

0 

1 

2 

3 

O  O'’  »  3e 

3 

-46 

12 

4.6670 

0 

1 

3 

1 

0  -  O*'  ♦  4e 

4 

-55 

20 

7.1872 

0 

1 

5 

O  -  O''’  ♦  5e 

6 

-62 

30 

10.1113 

0 

1 

5 

6 

0  -  O-'’  *  6c 

6 

-100 

42 

13.4128 

0 

1 

6 

7 

O  O*’  f  7e 

7 

-100 

56 

17.0724 

0 

1 

7 

8 

0  O'*  '  8c 

t> 

-100 

72 

21.0724 

0 

1 

8 

9 

H  H'  .  e 

1 

-20 

1 

2 

1 

0 

1 

10 

20  0; 

1 

100 

0 

0 

0 

2 

0 

11 

20  Of  '  e 

0 

55 

2 

1 

0 

2 

1 

12 

2H  lb 

-1 

20 

0 

0 

2 

0 

0 

iS 

?H  H;  ’  «• 

4 

2 

1 

•> 

0 

1 

1-1 

2II  •  O  -  H;0 

-2 

145 

0 

0 

2 

1 

0 

IS 

0  -  H  OH 

-1 

80 

0 

0 

1 

1 

0 

16 

0  -  H e  OH" 

-2 

no 

0 

0 

1 

1 

-1 

17 

0  +  H  OK'  ♦  0 

0 

50 

2 

1 

0 

1 

1 

13 

20  ’  2H  =  HjOj 

-3 

0 

0 

2 

2 

0 

19 

O  +  e  -  O" 

-1 

15 

0 

0 

-1 

0 

-1 

20 

H  ♦  e  H" 

-1 

5 

u 

c 

“  A 

1 

-1 

21 

20  =  o;  2e 

-1 

40 

6 

2.58V6 

1 

0 

2 

It  will  be  noted  that  the  reaction  equations  are  arranged  in  such  a  way  that  the  first  constituent 
on  the  right-hand  side  is  that  one  for  which  the  concentration  n;  is  given  by  Eq.  1.27.  (For  this 
constituent,  1  is  equal  to  j.)  There  are  no  such  equations  (or  O,  H,  and  e.  since  these  are  de¬ 
termined  by  Eqs.  1.24  and  1.26.  In  all,  we  have  21  equations  (or  the  concentrations  of  24 
compounds. 

Solving  these  equations  numerically  Is  extremely  difficult  because  the  equilibrium  constant, 
Eq.  1.23,  depends  m  k  .ind  n^  (as  defined  in  Eq.  1.17)  which  can  be  enumerated  only  if  the  ni’s 
are  known.  The  problem  was  coded  for  a  fast  electronic  computer  using  an  Iteration  process, 
but.  owing  to  operational  difficulties,  only  a  small  part  of  the  computations  planned  could  be 
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r()iiip!i'l('‘l.  !i;  partioular,  P  was  not  possible  to  obtain  results  for  very  high  temperatures,  and 
the  rorrertioii  for  the  electrostatic  forces,  <t>  in  Eq.  1.23,  was  not  applied. 

Ncvpi'theles.3,  It  was  possible  to  obtain  results  in  the  most  important  region  as  illustrated 
in  Fii;.  1.1.  The  machine  calculations  gave  J,  njcJ|V(Rp  Oj),  ng,  k,  and  S  -  Sj  directly  as 
functions  of  p,  v,  and  T.  y  was  computed  by  hand  from  Eq.  1.21  after  eraphical  differentiation 
of  p.g  .tnd  k.  The  shock  parameters  p„  v,,  and  i  j  =  yi(V|)  are  readily  obtained  by  graphical 
interpolation.  The  results  are  listed  in  Table  3.1. 


1.11  CALCULATIONS  USING  THE  THOMAS- FERMI-DIK AC  MODEL,  REGION  III 


For  temperatures  lover  than  those  considered  in  the  previous  section,  the  water 
molecule  remains  undlssoclatcd.  The  HKW  equation  of  state  is  still  valid  under  these  condi¬ 
tions,  out  an  ''ntirely  different  approach  is  possible  here,  namely,  the  use  of  a  statistical  model 
for  the  H2O  molecule. 

The  statistical  model  repla>  as  the  distinct  electron  orbitals  of  an  atom  by  a  continuous 
electron  cloud  which  ‘.s  treated  as  a  degenerate  Fermi  gas. 

This  model  Is  iiarticularly  suited  for  the  description  of  extremely  dense  matter,  and  it  has 
been  used  to  obtain  information  applicable  to  astrophyslcal  problems.  By  interpolations  be¬ 
tween  Bridgman's  experimental  data  and  the  results  of  theories  using  this  model,  it  was  also 
possible  to  obtain  information  of  interest  in  geophysics. 

The  same  approach  was  made  by  J.  H.  Rosenbaum  for  water  in  a  calculation  made  spe¬ 
cifically  for  this  project.  The  following  highly  Ideallxed  structure  was  assumed  for  water:  Each 
oxygen  molecule  is  surrounded  symmetrically  by  a  certain  number  of  hydrogen  molecules,  M. 
Two  of  these  belong  to  the  molecule  considered;  the  rest  belong  tc  the  nearest  neighbors.  The 
me'Uum  is  thus  represented  by  means  of  polyhedrons,  each  having  an  oxygen  ’.ucleus  at  Its 
center  anc!  M  hydrogen  nuclei  on  Us  surface.  Each  polyhedron  is  approximated  by  a  sphere,  and 
the  hydrogen  nuclei  are  smeared  over  the  surface  of  this  sphere.  To  be  electrically  neutral, 
such  a  sphere  must  contain  10  electrons,  eight  corresponding  to  the  nuclear  charge  of  the 
oxygen  atom  and  two  to  that  of  the  two  hydrogen  atoms.  If  the  radius  of  this  sphere  is  A,  the 
specific  volume  Is 


-  l!  (fl’  i!!.4v . 

m 


(1.28) 


where  N,\v  Avogadro’s  number  and  m  -  molecular  weight.  The  pressure  was  found  from  the 
virial  theorem  to  lie: 


P 


-  (1  -  c) 


(1.29} 


where  e 
n 

X 

l<0 

Z 

N 

c 


charge  of  the  electron 

Thomas- Fermi  unit 

reduced  radiu.s  -  |i<tt 

exchange  correction  -  ,to(Z) 

charge  nl  the  oxygen  nucleus  at  the  center  =  8 

charge  of  the  hydroge.i  nuclei  on  the  surface  of  the  sphere  =  2 

a  slowly  changing  function  of  M  (see  Table  I  in  NAVORO  Report  3847) 


The  potential,  'V,  of  the  electron  gas  Is  determined  by  the  TFD  equation 


with  the  following  Imundary  conditions 


(1.30) 
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Fig.  1.1 — lUnWne-HugorUot  leii'ioiu  fex-  w»ter. 


'!  (0)  -  1 


't'(x)  =  x+'Cx)  -  —  at  the  surface  of  the  sphere 

The  primes  Indicate  differentiation  with  respect  to  x.  I’ht  first  boundary  condition  yields  the 
potential  dun  to  ‘he  oxygen  nucleus  at  the  center;  the.  second  condition  ensures  that  H  positive 
charges  are  distributed  over  the  surface  of  the  sphere  and  that  the  whole  sphere  is  electrically 
neutral,  l.e.,  has  Z  N  electrons  inside  the  sphere.  A  more  complete  discussion  is  found  in 
NAVORD  Report  3847,  where  an  exliaustive  ll'jt  of  references  is  given.  The  results  of  these 
calculrtlnns  can  be  summarized  as  follows; 


V, 

P,  atm 

P,  atm 

cc/g 

(for  M  =  8) 

(lor  M  12) 

0.0988 

4.52  X  10’ 

4.55  X  10’ 

0.153 

1.78x10’ 

1.80  V  10’ 

0.269 

5.04  X I0‘ 

5.15  X  10* 

0.700 

4.44  V  10’ 

4.73  xic’ 

It  appears  that  the  pressure  is  relatively  insensitive  to  the  choice  of  M,  which  depends  on  the 
number  of  nearest  neighbors.  This  is  fortunate  because  the  actual  structure  of  water  at  the 
high  temperatures  and,  therefore.  M  are  unknown. 

The  pressures  given  above  refer  to  absolute  zero  temperature.  For  higher  temperatures 
the  “thermal"  pressure  Px  must  be  added.  In  NAVORD  Report  3847  the  following  e.xprcseion  is 
derived: 


This  can  be  evaluated  easily  lor  the  degenerate  electron  gas  and  leads  to  the  simple  formula 

p.r  =  2r.35  2  (1.3S) 

This  equation  of  state  (namely,  Eqs.  1.28,  1.29,  and  1.33  combined)  is  valid  only  for  extremely 
high  pressu"es.  Data  »•>'  Region  ID  can  be  obtained  by  interpolations  between  these  extremely 
high  pressures  and  Region  IV,  for  which  experimental  data  exist.  Figure  1.2  shows  such  an 
interpolation.  It  is  now  possible  to  derive  the  shock-front  parameters  as  well  as  the  other  data 
by  the  methods  described  In  the  previous  paragraphs.  These  data  apply  only  for  those  tempera¬ 
tures  where  no  substantial  dissociation  of  the  water  molecule  takes  place.  The  final  results  are 
shown  in  Fig.  1.1;  pi  and  n  as  obtained  from  the  TFD  model  are  lower  by  a  small  amount  than 
those  obtained  from  the  calculations  with  the  HKW  equation.  But,  in  general,  the  agreement  is 
as  good  as  can  be  expected,  considering  the  approximate  nature  oi  ootn  approaches.  Of  par¬ 
ticular  interest  is  the  agreement  in  the  temperature;  because,  according  to  a  theory  of  Jones,* 
even  a  crude  evaluation  of  detonation  rates  Is  expected  to  give  results  of  fair  accuracy  for  the 
p-v  data,  whereas  good  results  for  the  temperature  can  be  obtained  only  If  the  form  of  the 
equation  of  state  is  correct.  Our  results  may,  with  a  certain  reservation,  Indicate  that  the 
form  of  the  HKW  equation  is  approximately  correct. 

1.7  INTERMEDIATE-PRESSURE  RANGE,  REGION  IV 

The  equation  of  state  at  pressures  of  about  100  kb  and  below  has  been  treated  in  various 
publications  (see.  for  Instance,  references  6  to  8).  We  have  used  the  isentropic  Talt  equation 
as  proposed  by  Kirkwood  and  Bethe,  Eq.  2.2  of  reference  0,  to  map  out  y  and  J  in  the  region 
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Fig.  1.2 — CalcuUtloni  uring  the  TFD  mo&.:l. 
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between  the  shock  front  and  the  saturation  line.  It  was  necessary  to  use  variable  coefficients 
and  exponents  in  this  equation.  This  Talt  equation  failed  completely  in  and  near  the  region 
where  the  HKW  equation  was  used.  The  results  are  shown  In  Figs.  1.3  and  1.4. 

1.8  RANGE  OF  VERY  LOW  PRESSURES,  REGION  V 

For  the  sake  of  completeness,  the  relations  which  hold  for  low  pressures,  where  the  shock 
wave  behaves  almost  like  an  acoustic  wave,  are  given  below.  If  the  sound  velocity  c  changes 
linearly  with  pressure,  i.e.. 


C  =  Cj(l  +  Jp), 

(1.34) 

the  following  relations  can  be  derived: 

Vj  -  V  _  p  1 

V  Pocl  1  *  tp 

(1.35) 

(1.38) 

(1.37) 

For  very  low  pressures  the  following  first-order  approximations  can  be  derived  for  the  shock 
front: 

(1.38) 

(1.39) 

where 

a,  =  fp,c| 

a,  =  1 

(1.40) 

Using  these  first-order  approximations,  we  obtain  for  the  dissipated  enthalpy  increment  h  (see 
definition  in  Eq.  1.14) 

llmp,h=.jP.(P‘p/'). 

(1.41) 

In  these  equations  we  have  used  the  relation  Vj  =  l/pj. 

Strictly  speaking,  the  magnitude  i  as  defined  by  Eq.  1.34  refers  to  the  change  of  sound 
velocity  along  an  Isentropic.  For  low  pressures,  however.  It  is  permissible  to  use  values  of  C 
which  are  obtained  by  measuring  sound  velocities  as  functions  of  pressure  at  constant  tempera¬ 
ture.  Experiments  recently  completed  at  the  Naval  Ordnance  Laboratory"  give 

«  =  0.106  kb-‘. 

This  value  holds  for  sea  water  as  well  as  for  fresh  water  over  a  wide  temperature  range  (0  to 
40*0  which  includes  the  range  tn  which  we  are  interested.  For  a,  we  have  used 
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ill  ^  2.J7fi. 


Thin  value  holds  (or  fresh  water  at  20*C  and  also  (or  sea  water  at  about  8*C.  The  Utter  ciee 
corresponds  lo  the  conditions  at  Wigwam,  whereas  the  former  ease  U  that  for  which  the 
thermodynamic  calculations  were  made.  It  seems  reasonable  to  asstiiie  that,  for  high  pres- 
surcB  where  the  pertinent  Information  Is  lacking,  our  calculations  which  were  made  for  fresh 
water  ut  20°C  hold  adequately  for  sea  water  at  a  temperature  of  8*C.  Since  In  the  shock>wave 
calculations  only  the  ratio  Pi/pt  occurs  and  the  absolute  value  of  the  density  of  the  ambient 
medium  does  not  appear,  these  calculations  are  relatively  Insensitive  to  the  ambient  tempera* 
ture. 


1.9  RESULTS 

The  (Inal  results  (or  the  shock-front  parameters  are  listed  In  Table  3,1  and  are  shown  in 
rig.  1.1.  The  Isentroplc  exponent  y  and  the  reduced  internal  energy  J  for  the  whole  region  of 
our  Interest  are  presented  In  graphical  form  in  Pigs.  1.3  and  1.4. 

In  general,  thermodynamic  (unctions  obtained  from  equilibrium  calculations  are  not  smooth. 
With  Increasing  temperature,  different  dissociation  and  ionisation  levels  are  reached,  and 
therefore  the  number  of  particles  changes  discontlnuously;  this  resulU  in  irregularities  In  the 
graphical  representation  of  such  functions.  For  practical  reasons  we  have  used  smoothed  daU 
In  the  hydrodynamic  calculations,  and  Table  3.1  gives  such  data  with  these  Irregularttlrs  elimi¬ 
nated.  In  the  region  where  there  la  an  overlap  in  the  data  given  by  the  TFD  model  and  the  HKW 
equation,  we  have  given  the  latter  full  preference  (except  for  the  isentroplc  exponent  >),  al¬ 
though  the  differences  are  not  very  large.  The  HKW  p-v  data  are  in  excellent  agreement  with 
the  experimental  results  of  Walsh  and  Rice,**  who  measured  propagation  and  particle  velocities 
of  strong  shock  waves  in  water.  Their  results  are  indicated  in  Fig.  l.i. 

The  results  of  our  equilibrium  calculations  of  the  Isentroplc  exponent  y  show  a  greater 
scatter  than  that  obtained  (or  the  other  thermodynamic  variables.  This  scatter  in  y  may  be 
attributed  to  variations  caused  by  the  Inclusion  of  excited  states  and  to  uncertainties  Introduced 
by  the  graphical  differentiations  required  (see  Eq.  1.31).  For  y  to  vary  smoothly  from  the  low- 
pressure  region  to  the  value  1.666  at  infinite  pressure,  it  was  necessary  to  use  an  average  be¬ 
tween  the  data  obtained  from  the  TFD  model  anu  the  HKW  valuea  (see  Fig.  1.1). 

y  varies  from  the  value  1.666  (at  Infinite  pressure)  to  infinity  (In  the  limit  of  aero  excess 
pressure).  Attempts  were  made  to  Introduce  new  variables  so  that  a  more  well-behaved  func¬ 
tion  could  be  used  In  place  of  y.  It  was  not  possible  to  obtain  substantially  simplified  relations, 
and  thus  11  was  r  cessary  to  use  the  calculated  values  of  y  in  tabular  or  graphical  form. 


26 


SECRET 


RESTRICTED  DATA 


CHAPTER  2 


CALCULATION  OF  SHOCK-WAVE  PHENOMENA 


2.1  INTRODUCTION 

The  calculations  of  the  shock-wave  phenomena  produced  by  an  atomic  buret  were  based  on 
the  Idealtzlng  concept  of  a  point  explosion,  i.e.,  an  Infinitely  small  explosive  charge  delivers  a 
finite  amount  of  energy  Instantaneously.  At  the  first  moment,  such  a  concentrated  energy  dis¬ 
charge  produces  Infinitely  high  pressures  and  temperatures.  An  Instant  later,  these  become 
finite  and  subsequently  run  through  the  range  which  would  occur  In  an  atomic  explosion.  It  Is 
this  later  range  which  wu  shall  try  to  calculate  here,  since  Initially  the  model  of  a  point  ex¬ 
plosion  is  unrealistic  even  for  atomic  explosions,  which  produce  high  but  finite  temperatures 
and  do  not  release  their  energy  at  a  single  point. 

An  impoi  iant  consideration  lor  atomic  explosions  Is  the  question  of  the  energy  transport  by 
raoiation.  A  study  of  this  question  by  Weber'‘  showed  that,  for  underwater  explosions,  radiation 
phenomena  are  much  less  Important  than  for  air  bursts.  In  our  calculations  we  have  neglected 
radiation  entirely  and  have  assumed  that  the  energy  of  the  burst  Is  transmitted  from  one  par¬ 
ticle  to  the  next  by  pressure  forces  only.  This  Is  a  satisfactory  assumption  once  the  shock 
Irunt  has  traveled  a  certain  distance,  although  It  Is  a  pcor  assumption  for  the  early  stage  of  the 
explosion  and  lor  the  region  near  the  center  for  later  stages.  Since  we  are  not  concerned  with 
these  phenomena  here  and  since  they  do  not  affect  the  results  of  Interest  to  us,  the  neglect  of 
radiation  seems  appropriate. 

The  calcuiatlons  described  in  this  chapter  .-leal  with  the  task  of  Integrating  the  spherical 
blast-wave  equations  with  variable  Isentroplc  exponent  y.  They  had  to  be  carried  out  without  Ute 
benefit  of  an  electronic  computer,  and  so  certain  approoclmatlcns  were  necessary.  However, 
none  of  the  approximations  made  seem  to  be  of  a  serious  nature.  The  greatest  advantages  In 
tractablUty  were  achieved  by  the  introduction  of  reduced  variables  and  a  corresponding  trans¬ 
formation  of  the  hydrodynamic  equations. 


2.2  ENERGY  EQUATION  AND  PRESSURE-DISTANCE  RELATION 


The  total  energy  contained  within  the  sphere  bounded  by  the  shock  front  is  constant  and  is 
equal  to  Q,  the  hydrodynamic  yield  of  the  weapon.  All  particles  within  this  sphere  have  kinetic 
as  well  as  Internal  energy.  The  law  of  conservation  of  energy  may  be  written: 


Q-  '*-'  r  fe’+P(E-E,) 
•'«  I  ^ 


r*  dr  =  constant , 


(2.1) 
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whore  E  denotes  >he  internal  energy  per  unit  mass  for  partirles  inside  the  sphere  of  radius  r| 
and  L'd  denotes  the  Internal  energy  per  unit  mass  of  the  ambient  water.  We  transform  this 
eauatlon  by  introduction  of  the  following  reduced  variables: 


tl  ■-  reduced  pressure  = 


pressure _ _  p 

shocl^wave  peak  pressure  Pt 

velocity  _  u 


9>  =  reduced  particle  velocity  •-  —  ,  ^ 

•'  peak  velocity  U| 


X  =  reduced  density  = 


4  =  reduced  distance  = 


density  _  p 
peak  density  p, 

radial  distance 


radial  distance  of  shock  front  r| 


E  -  Et 

J  =  reduced  Internal  energy  = - -  p  (see  Eq.  1.4). 


With  these  new  variables  and  with  the  use  of  the  Ranklne-Hugonlot  conditions  given  in  Eqs.  1.1 
and  1.3,  Eq.  2.1  takes  the  form 


/N  _  Pi  -  Pj  _ 

where  ih  is  the  reduced  total  energy  defined  by 


(2.2) 


(2.3) 


Since  Pi  is  a  function  of  p,,  Eq.  2.2  provides  a  convenient  relation  tor  the  peak  pressure  pt  as 
a  function  of  distance  r|  lor  any  given  yield  Q,  once  qi  is  known. 


2.3  HYDRODYNAMIC  EQUATIONS 

To  determine  r/|  we  have  to  find  solutions  of  the  partial  differential  equations  of  the 
spherical  fluid  motion 


u,  +  uur  p  Pr  =  ® 

(2.4) 

P|  +  uPr  +  P“r  =  0 

(2.5) 

P|  +  »‘Pr  =  'HP>P)(Pt+''Pr) 

(2.0) 

and  their  boundary  conditions,  which  are  the  Rankine-Hugoniot  conditions,  Eqs.  1.1  to  1.3.  In 
the  above  equations,  p  denotes  the  excess  pressure  above  the  hydrostatic  pressure  P^.  y  is  the 
Isentroplc  exponent  discussed  in  Sec.  1.2. 

By  introduction  of  tli-;  reduced  variables  p,  x>  and  and  the  Ranklne-Hugonlot  conditions, 
the  hydrodynamic  equations,  Eqs.  2.4  to  2.8,  can  be  brought  into  the  following  form; 

J_  Po  +  L  +  G  * 

^(vxPt-Pt  ’  *  a  In  Pi  ,2 

'T  Pi  <P  '  ■ 

Pt  -  P«  “  £ 
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(2.8) 


4.r  *  X 

£X  _  X  4  *  *  8  In  Pi 

»{"«  Pi  y 

Pi  “  Pt  I 

/3«)  .  2«>\  8  ln» 

8»  _  »  ^  j  a  In  Pi 

a{  7  Pi  » 

Pi-Pa"  4 


_3  /_Pi_\ 

H.  -  H.(p..«  =  ^ i|Pt - - 

yR„pi-P.' 


G|  =  G,(p,./))=^i^  =  :^H, 
u,  dt  yRj, 

L.  =  L,(p,.«  =iii^= »!  A  +-1-  -et_\ 

u,  dt  2\  yjufPi-PeJ 


r  =  r(Pi,x,W  yRH  “  yRH  (Pi> 

j_Pi-P<rilni;i 
p,  'dlnp,' 


It  wUl  be  noted  that  the  Ume  variable  haa  been  ollmlnated  from  these  equations  and  ma]i  be 
regarded  as  replaced  by  th«  peak  density  p,.  The  boundary  conditions  are  siupiy 


S’  - 1) 

(P  =  0) 

X  «  =  1 

^-0 

V-  =  i! 

(The  latter  two  boundary  conditions  are  obtained  from  symmetry  considerations  at  the  center  of 
the  sphere  and  from  the  hydrody.iamic  Eq.  2.4.)  It  follows  that  the  derivatives  of  x,  and  ^ 
with  respect  to  P|  vanish  for  {  =  1.  Thus,  for  {  -  1,  Eqs.  2.7  to  2.9  become 

/a^X  ^2yiJj9^ 

\a?/(-i  _£i__y, 

Pi-Pi 

X\  =  {%)  =ei^‘(?i  +  2*G,)  (2.1?) 

W/^.l  Pc 

=  (If)  2)  +  (*-18) 

These  differential  quotients  are  functions  of  P|  and  fi,  where  d  Is  an  as  yet  undetermined  func* 
tion  of  Pi,  as  indicated  by  Eq.  2.14. 
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2,4  NATURE  OF  THE  SOLUTION 


In  ordei  to  see  how  J  Is  related  to  the  solution  of  our  problem,  we  use  the  followtng  con¬ 
sideration:  If  r(Pip)  Is  an  analytic  function  at  the  shock  front,  we  can  calculate  the  higher 

derivatives  </>'{,  tpf' . xWxTi  •  •  •<  VTi  tif  •  •  •  •  •  All  are  functions  of  p|  and  p.  Assuming  <p, 

X,  and  i  to  be  analytic,  we  can  make  a  Taylcr  expansion  around  {  -  1: 


(nU.Pi.P)  =  1  -  -  0  ^  (1  -  0*-  ^  (!-«)»  +  ... 

(2.19) 

>M,PuP)  =  1  -  xid  -  0 + ^  a  -  f)*-  (1  -  i)*  + . .  . 

(2.20) 

4(f,p„(i)  =  (i-4)‘+ ... 

(2.21) 

Introduction  of  these  expansions  into  the  integral,  Eq.  2.3,  tor  the  reduced  total  energy  it) 
yields  >;i  as  a  function  of  p|  and  0.  This,  by  virtue  of  Eq,  2.14,  is  an  ordinary  differential  equa¬ 
tion  of  the  first  order  for  iji.  Therefore  this  approach  has  reduced  the  three  partial  differential 
equations  to  one  ordinary  differential  equation. 

7cr  practlc.'tl  calculations  this  method  is  not  suitable  because  the  coefficients  of  the  hlgh- 
crder  terms  in  the  expansions,  Fqs.  2.12  to  2.21,  are  too  complicated.  Therefore  a  combination 
of  expansions  around  f  =  0  and  {  =  1  was  used  in  this  project,  as  will  be  described  In  the  follow¬ 
ing  sections.  The  differential  equation  for  i})  was  solved  graphically  by  the  method  of  Isoclines. 

2.S  TAYLOR  CASE 

In  Sec.  1.3  we  have  seen  tliat,  for  Infinitely  high  pressures,  the  density  pi  directly  behind 
the  shock  front  approaches  a  finite  value,  pi  =  4p|.  We  also  have  y  %  and  l/ygH  =  0.  In  this 
case  the  three  partial  differential  equations,  Eqs.  2,7  to  2.2,  become  ordinary  differential 
equations.  This  Is  the  well  known  Taylor**  solution  for  strong  point  blast  waves. 

We  shall  assume  that  the  Taylor  solution  gives  the  density,  velocity,  and  pressure  distri¬ 
bution  within  the  Infinitesimal  sphere  at  t  -  0.  This  Is  an  Idealising  approximation  because 
actual  explosions,  even  from  concentrated  energy  sources,  canno,  ,ave  such  distributiunif  from 
the  beginning.  However,  It  seems  safe  to  assume  that  the  actual  distributions  converge  quickly 
to  those  which  are  calculated  with  the  assumption  of  an  Initial  Taylor  solution. 

Since  Infinitely  high  temperatures  are  assumed  to  prevail  throughout  the  sphere  of  dis¬ 
turbance  at  t  =  0,  the  isentroplc  exponent  y  is  constant  within  this  sphere  and  equals  '/|.  For 
constant  y  the  differential  equations  can  be  solved  In  closed  form  (see  NAVORD  Report  4182). 
The  solution  Is  shown  In  Fig.  2.1  In  the  curves  labeled  pt  =•  4pi.  It  Is  Interesting  that  the 
density  vanishes  at  the  center  of  the  explosion.  The  same  bolds  true  at  later  times  when  the 
shocK  front  has  traveled  to  larger  distances. 

2.8  EXI‘.*J»SIONS  AROUND  (  =  0 


The  Taylor  solution  can  be  given  not  only  in  closed  form  (which  Is  discussed  in  NAVORD 
Report  4182)  but  also  in  the  form  of  Infinite  series; 

^  „  .m**  «  .2m+5 

X  =  Bj?"*  +  +  (2.22) 
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whurt 


4.5 


Explicit  exprcs.sloiis  for  the  first  few  coefficients  are  also  given  in  NAVORD  Report  4182.  The 
numcricul  va'ues  nre  =  0.306,  Aj  =  0.320,  Bj  =  0.343,  C,  =  C.800,  and  Cj  =  0.229. 

In  NAVOno  Report  4183  it  is  shown  that  such  expansions  may  be  used  to  express  the 
desired  soiutions  of  the  partini  differential  etjuations,  Eqs.  2.7  to  2.14.  For  this  general  case, 
that  is,  i/viiH  *  0,  the  coefficients  in  Eq.  2.22  are  functions  of  pj.  Substituting  these  expansions 
into  the  partial  differential  equations  and  equating  like  powers  of  {  yields  a  set  of  ordinary 
dilferential  equations  in  the  coefficients  Ai,  Bj,  and  Ci.  The  first  three  of  these  are 

3nCo  t  H,  t  O,  =  0  (2.23) 


ni 


^Pl-Po 


/ 


SCj  —  Gf  —  G| 


d  In  B, 
d  In  Pi 


=  0 


(2.24) 


Po  ni  t-  2 
PI  -  PO  I'oCo 


A,  + 


L, - +  Cl,  +  G. 

P|-Po 


d  In  Cq 
dlr.p, 


0. 


(2.25) 


These  are  three  equations  In  the  four  unknowns  A^,  A|,  Bt,  and  C^.  Foi  each  additional  differen¬ 
tial  equation  added  to  this  system,  one  more  unknown  is  introduced.  Therefore  the  expansion 
around  f  -  0  is  undetermined  unless  additional  information  is  incorporated. 

As  described  in  NAVORD  Report  4184,  this  may  be  accomplished  in  several  different 
ways.  The  most  successful  method  tried  uses  Eq.  2.7.  After  solving  for  and  integrating 
from  .e  0  to  £  1.  wo  obtain 


Ao 


.  I-'.'',  f 

Po 


Pi  -  Po 
Po 


fl  /b^ 

Jo  W 


+  L|  +  G| 


9  In  y 
0  In  Pi 


(2.26) 


2./  POLYNOMIALS 

.Since  the  solutions  can  be  expressed  by  means  of  expansions  around  £  =  1  (Eqs.  2.19  to 
2.21)  and  around  £  -  0  (Eq.  2.22),  one  may  construct  approximate  solutions  by  using  relatively 
few  terms  of  each  and  merging  these  together  in  the  intermediate  range.  One  of  the  simplest 
way.s  to  do  this  is  to  use  ihc  expansion  a.'ound  £  =  0  for  the  entire  range  of  £,  but  to  determine 
a  curtain  number  of  tbc  coefficients  in  such  a  way  that  the  boundary  conditions,  Eq.  2.15,  are 
fulfilled  as  well  as  the  behavior  near  £  =  1,  as  given  by  the  first-order  derivatives,  Eqs.  2.16 
to  2.18,  and  any  higher  order  derivatives  which  one  may  wish  to  Include. 

If  only  first-order  terms  are  used,  we  employ  a  four-term  polynomial  for  x  and  a  three- 
term  polynomial  for  and 

V  Cjf  ^  C|£  I 

\  -.Bofin  +  B|£2'n*Z+ B2£®"'*‘‘+ B|£^"’**.  (2,27) 

£  --  Ao  +  Aif'"’^  +  a,£2'"*4. 

There  arc  then  iO  coefficients  to  be  determined,  for  which  we  have  the  following  conditions: 

1.  Three  boundary  conditions,  Eq.  2.15. 

2.  Three  first-order  derivatives  at  £  =  1,  Eqs.  2.16  to  2.18. 

3.  IWo  first-order  differential  equations,  Eqs.  2.23  and  2.24. 
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One  further  condition  can  be  obtained  from  the  fact  that  the  mass  of  the  sphere  which  Is 
bound  by  the  shock  front  Is  equal  to  that  of  an  equally  large  sphere  containing  water  at  normal 
density: 


♦"Pir!  jT*  xl*  «1«  =  y  P»ri 


B|  ^  B|  ^  B;  ^  Bj  _  _Po_ 
m  +  3  2m  +  5  3m  +  7  4m  +  9  3pj 

The  tv'nth  condition  is  the  relation  given  as  Eq.  2.26.  Introducing  the  polynomials  in  Eq.  2.27 
Into  Eq.  2.26  yields,  with  the  use  of  the  abbreviation 


6.5bj=^j^  J  =  l,  2.  3 .  (2.29) 

the  following  relation  for  A«: 

Ao  =  1  +  L,|c,(b,  -  2bi  +  b,)  +  (14  -  ft)  ^2^  +  b,j 

-  ^  C,(b,  -  2bi  +  b,)  +  (14  -  ipi)  +b, j 

+  £l^  [C|b,  +  8.5CoC,b,  +  (7.5C5  +  15C,C,)b,  +  21.5C,Cjb4  +  MClbs) 

+  CLZ-CJ  G,  (b,  -  2bj  +  b.) - ^  f  .  (2.30) 

p,  ''‘l«'«nPt  '  '  ’  dlnpil  6,5  ;J 

This  Is  an  ordinary  dl/ferential  equation  for  C(  as  a  (unction  of  pj.  Thus  we  have  three  simul¬ 
taneous  ordinary  differential  equations  which  determine  A«,  I3g,  and  Cg,  namely,  Eqs.  2.23, 

2.24,  and  2.30.  The  other  coefficients  are  given  as  follows: 

A,  =  ll-2A,-tf',/6.5 
Ag  =  1  -  Ag  -  /•  t 

B,  =  -  5,6B,  +  — -  7,3728  +  0.1657x', 

B,  =  8.2Bc  -  ^  18.4379  -  0.4852xi  (2-31) 

*  Pi 

B,  =  -  3.6B,  +  _  10.0651  +  0.3165x', 

C,  =  ^-2C,-«.',/6.5 


C,  =  1-C,-C,. 
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The  arcuracy  of  these  polynomials  has  been  tested  for  the  Taylor  case,  where  the  exact 
solullon  is  known.  There  was  agreement  up  to  the  third  decimal  place,  which  showed  that  It  Is 
not  necessary  to  include  any  more  terms  In  the  polynomials.  This  seemed  particularly  Justified 
in  view  of  the  fact  that  v.’c  are  not  sc  much  Interested  in  the  reduced  functions  <p,  y,  and  if  as  In 
the  Integral  formed  from  them.  Small  Inaccuracies  or  approximations  used  under  such  an 
integral  usually  have  only  a  minor  Influence  on  the  accuracy  of  the  final  result. 

2.8  DETERMINATION  OF  p 

In  all  relaticns  derived  so  far,  p,  defined  by  Eq.  2.14,  occura  as  an  unknown  function  of 
Pi.  Only  those  «>’s.  x’s,  and  tl*s  which  involve  a  p  satisfying  equation,  Eq.  2.14,  are  solutions  of 
our  problem. 

The  method  used  in  this  project  determines  P  by  means  of  an  Interpolation.  First  aome 
arbitrary  and  constant  values  are  assumed  for  p.  Then  (p,  y,  and  are  determined  by  the 
methods  described  above.  Subsequently,  i)i  Is  calculated  using  Eq.  2.3.  When  In  rii  is  plotted  vs 
In  (pf  -  pt)/pc,  the  inclination  of  the  r/,  curve  we  seek  must  have  the  value  of  P  wMch  waa  uaed 
to  calculate  this  curve.  Such  a  curve  is  readily  found  in  the  tame  way  as  the  graphlcsl  solution 
of  differential  equations  by  the  method  of  isoclines.  Figure  3.2  illustrates  the  procedure. 

Since  It  turned  out  that  r/t  was  not  sensitive  to  the  value  assumed  for  p,  an  interesting  con¬ 
clusion  may  be  made  on  the  nature  of  the  solution.  If  were  assumed  to  be  entirely  Independent 
uf  the  assumed  value  of  p,  Eq.  2.14  would  not  be  needed  at  all  to  arrive  at  a  solution,  and  this 
would  eliminate  the  arbitrary  integration  constant  of  this  differential  equation,  l.e.,  we  would 
have  a  singular  solution,  Independent  of  the  Initial  condition.  It  appears  that  this  situation  is 
approximately  true  for  the  solutions  we  have  obtained.  It  therefore  seems  possible  to  start  the 
calculation  with  different  Initial  Taylor  distributions,  l.e.,  with  different  values  of  y’,  such  as 
1.4,  1.66,  or  7,  and  still  obtain  almcat  the  same  result  for  the  shock-wave  parameters  at  dis¬ 
tances  where  pressure  measurements  can  be  made. 


2.9  LOW-PRESSURE  RANGE 

When  the  shock  front  has  propagated  to  large  distances  and  the  pressure  Is  low,  the  wave 
behaves  nearly  like  an  acoustic  wave.  The  peak  pressure  of  such  waves  decreases  approxi¬ 
mately  proportionally  to  the  dhstance.  For  Instance,  Kirkwood  and  Bethc*  have  found  the 
following  asyinplulic  behavior  uf  the  weak  shock  waves; 


r,  (in  r,/r,)'j’ 

where  c  Is  a  constant  and  rg  is  the  reference  radius.  From  this  we  obtain 


(2.32) 


d  Inpi  _  j  ^  1 _ 1 

~  d  In  Tj '  2  In  r/rg 

This  magnitude  approaches  unity  slowly  as  the  shock-front  distance  increases  to  Infinity. 
Going  back  to  Eq.  2.10,  we  find 


(2.38) 


d  In 


£i  =  Hh.: 


dlnr,  u,-'  <  .  J_  _£i_  +  0) 

Vrh  Pi  “  P* 

or  for  low  pressure,  with  the  use  of  the  approximation,  Eq.  1.38, 


a  in  Tj  J 


(2.34) 


(2.35) 
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ConipariHun  o{  Eqs.  2.:*S  and  2.33  shows  that,  for  low  pressure,  /I  must  slowly  .approach  unity. 
From  the  definition  of  fq.  2.14,  it  can  be  seen  that  jji  must  tend  to  aero.  Thlf  gives  us 
qualitative  information  on  </)  ami  vi-  they  must  tend  to  vanish  for  most  values  oi  i  for  very  low 
shock  pressures,  or  else  the  integral  for  ;7|,  Eq.  2,3,  could  not  become  zero,  x  does  not  vanish, 
as  can  be  seen  from  the  average  density  condition,  Eq.  2.28;  in  fact,  x  approaches  unity  when 
ft  “  Pf 
Since 


</>  ~  1 
c  -  1 

a. id 

y  =  0  for  i  -  0, 

we  obtain  in  the  !lmi;  (or  zero  .shock  pressure: 


lor  t  --  1 


lim  <n  =  Urn  -  0  0  s  {  <  1 

**1  **>4 

llin  X  =  1  0  <  {  i-  1. 

Pi-rt 

There  is  a  discontinuous  drop  of  o  and  ^  from  1  to  0  at  {  =1  and  also  for  x  at  {  =  0. 

This  implies  that  the  derivatives  <o',  and  if'i  Increase  to  Infinity  when  Pt  —  pt,  but  y', 
vanls.hes.  Indceu,  Eqs.  2.16  to  2.18  become,  with  the  use  of  Eqs.  1.38  and  1.39, 


lim  =  Urn  t\ 

ft-Pt  pi-pt 


2p,  I  -  If 
Pi-Pt  aj2  +  ()’ 


lim  y; 


i  *  -■< 

Tj  2  +  d 


0. 


(2.36) 


Thus  (!  -  d)/(pt  -  Pi)  approaches  •>  as  p,  approaches  pg.  It  is  difficult  to  evaluate  the  integral  Vi 
for  these  conditions,  because  J  Increases  rapidly  with  decreasing  (  for  low  values  of  P|.  Basi' 
rally,  this  is  because  the  bubble  aiid  shock-wave  phenomena  become  spatially  separated.  Most 
of  the  shock-wave  e  icrgy  Is  concentrated  near  the  shock  front  and  can  be  easily  accounted  for 
by  the  methods  described  above.  But  this  Is  only  a  fraction  of  the  total  energy,  since  an  appre- 
cl.ible  amount  of  energy  Is  found  as  internal  energy  near  the  center.  This  latter  energy  pro¬ 
duces  the  phenomenc  of  the  pulsating  bubble  and  Is  not  conveniently  accounted  for  by  the  method 
described  in  Secs.  2.7  and  2.8. 

Therefore,  to  obtain  we  may  calculate  the  shock-wave  energy  and  try  to  find  a  relation 
between  total  energy  and  shock-wave  energy. 

At  the  shock  front,  energy  is  dissipated  because  of  the  Irreversible  shock  process  and  is 
converted  into  thermal  energy.  The  shock-wave  energy  Qg  Is  reduced  correspondingly.  Kirk¬ 
wood  and  Brinkley,  Eq.  4.29  in  reference  14.  have  formulated  this  phenomenon  as  follows: 


^  =  -  4rp,hr;,  (2.37) 

dr, 

where  Qg  =  4i:  (r(tll’  up  dt  and  h  dissipated  enthalpy  increment.  The  Integral  for  the 

shock-wave  energy  has  to  be  carried  along  the  path  of  the  fluid  pa.'ticle  which  is  reached  by  the 
shock  front  at  t  -  tfr,)  and  r  =  r,.  The  upper  limit  corresponds  to  a  long  time  compared  with 
the  duration  of  the  shock  wave  but  must  be  chosen  in  such  a  way  that  secondary  pulses  are  not 
included, 
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In  NAVORD  Repurt  4182  the  shock-wave  energy  Qg  is  expressed  by  a  volume  integral. 
Defining  the  reduced  shock-wave  energy  tjs  the  same  way  as  the  reduced  total  energy,  the 
following  expressions  have  been  found 


and 


lim  iJs=3  jf  +  (^d( 

p,~f,  °  •'0 


2¥>'l  +  X't  +  3  2«5  +  3 


(2.J8) 


(2.391 

(2.39a) 


where  the  derivation  of  Eq.  2.39a  is  explained  in  NAVORD  Report  4182.  With  Eq.  2.36,  this 
becomes 


3 

lim  =-  5  a, 

ei-Pi  °  2 


Pi  -  Po  2  *  <1 

Po  1  -  U  ‘ 


(2.40) 


This  equation  is  of  little  help  yet,  since  it  relates  the  reduced  shock-wave  energy  with  /9,  the 
derivative  of  the  reduced  total  energy.  Another  expression  for  Vs  he  obtained  from  Eq. 
2.37.  Introducing  Eq.  2,38  yields  after  a  few  manipulations 


'lc  = 


2hpJ 


1  + 


Pi 


Pi-Pt 


(1  +« 


Pi(Pi-P«)  f 

r, 


_£l_  (0o-fi) 


.ind,  with  the  use  of  Eqs.  1.38,  1.39,  and  1.41, 


(2.41) 


where 


Urn  iJs 

(>i  'pt 


a,  Pi  -  p,  2  +  0  ^ 
3  Pt  0s  ~  0 


0£ 


d  In  tjg 


d  In 


Pi  ~  P» 
P« 


(2.42) 


Combination  of  Eqs.  2.40  and  2.42  gives  tiie  following  differential  equation  for  r/g  (the  limit  sign 
is  omitted): 


3a2  Pi  -  Po(4/3)  +0S 
2'  p,  l-Ps 


Us  snliition  is 


Vs  = 


3a|  p.  —  Pc  /7 
2  p,  \3 


In 


^ - 1„-M5_), 

Pt- Pf  P\-Ptl 


(2.43) 


(2.43a) 


where  C  is  the  integration  constant.  Combination  of  this  solution  with  Eq.  2.40  yields  the  de¬ 
sired  asymptotic  expression  for  0\ 


lim  0  - 


3 

i  *  n 


(2.44) 
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The  .'iicnt  »[  this  slnipie  expression  is  that  it  shows  how  fl  can  be  extended  into  the  low- 
pressuie  repion  where  the  calculation  of  i}|  becomes  inconvenient.  The  integration  constant  c 
is  determined  by  the  results  for  jl  obtained  In  the  high-pressure  region;  0  must  merge  smoothly 
from  one  region  into  the  other. 


2.10  PEAK  APPROXi;.lATiON 

Another,  more  elaborate,  treatment  of  the  low-pressure  range  is  done  ‘'y  the  use  of  the 
nnsii  upnroxlmatlon.  By  "peak  approximation"  we  reier  to  the  attempts  of  various  authors  to 
describe  the  shuck-wave  prupagaiton  in  terms  of  shock  parameters  only.  This  can  be  dune  by 
making  certain  idealizing  assumptions,  for  instance,  assuming  a  similarity  restraint  on  the 
wave  shape.  Kirkwood  and  Bethe  as  well  as  Kirkwood  and  Brinkley  have  presented  such  shock- 
wave  theories,  and  in  their  application  to  underwater  explosions  they  assumed  an  exponential 
wave  shape.  These  theories  have  been  successful,  and  they  would  be  quite  appropriate  for  our 
case.  We  have  used  the  more  recent  theory  by  Snay  and  Matthias.* 

The  Snay-Matthlas  theory  gives  an  Interrelation  between  peak  pressure,  time  constant,  and 
the  profile  of  the  wave  which  is  characterized  by  the  shape  factor  i  This  magnitude  is  defined 
as  follows: 


where  f  is  unity  fur  an  exponential  wave  shape  and  zero  for  a  triangular  wave  shape.  The  theory 
leads  to  the  following  dlfferentiai  equations:* 

^  ^  +  Pii«  =  0  (2.46) 

^  *  P««*  (1  +  fPii)  =  0  (2,47) 

where  x  =  reduced  distance  =  r|/r» 

r^  =  reference  length  (will  cancel  in  this  analysis) 
a  3  time  factor  =  r«/  Vc, 

6  e  time  constant  =  -  p/(*p,'3t) 

C|  =  sound  velocity  at  p  =  0 

The  variable  coefficients  P|  of  these  differential  equations  are  obtained  from  the  three  hydro- 
dynamic  partial  differential  equations  and  are  functions  of  the  peak  pressure  P|  o.nly. 

Introducing  0  and  p|,  these  two  equations  can  be  transformed,  with  the  use  of  Eqs.  2.34  and 
1.12,  to: 


0  = 


W|  +  WjOX 


..  P*  -  PS  4 
-trh  — -1 


(2.48) 


and 


d  In  ox  -1  /  < 

- = - j  w, - a  -  W|OX  -  Wifax 

Pi  -  m  w,  ♦  W]ax\  *  ox 

Pi 


d  In 


)■ 


(2.49) 


*See  reference  1,  Eqs.  1.6  and  2.7,  where  the  symbol  ot/o*  Instead  of  f  is  used  to  desig¬ 
nate  the  shape  (actor. 
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where  w, 


*  Pi  5l! 

^RH  Pi  “  P*  Pi 

^  Pi  Pit 

^RH  Pi  ~  P*  Pi 


Wj  =  1  -  Fi4 

««  =  Pu 
w»  =  Pii 
i*»  =  P||Pi« 


anti 


:lL. 

C(|9 


Tne  variable  coefdclerits  W|  are  listed  in  Table  2.1.  They  are  nuw  functions  of  p|. 


(2.50) 


Table  2.1  —  VARIADLE  COEFFICIENTS  OF  EQS.  2.48,  2.48.  AND  2.SS 


Pl/Po 

W| 

Wj 

W> 

W4 

Wt 

Wi 

1.60 

0.3S23 

0.0375 

1.311 

4.45 

-0.0187 

0.098 

1.50 

0.3917 

0.0468 

1.282 

3.32 

-0.0182 

0.114 

1.40 

0.4578 

0,0607 

1.240 

2.33 

-o.otse 

0.128 

1.30 

0.5182 

0.0756 

1.190 

1.50 

0.0000 

0.140 

1.20 

0.S874 

0.0877 

1.124 

0.810 

0.0300 

0.142 

1.15 

0.6434 

0.0931 

1.092 

0.533 

0.0750 

0.129 

1.10 

0.7326 

0.0921 

1.067 

0.312 

0.0813 

0.110 

1.075 

0.8005 

0.0820 

1.052 

0.220 

0.0770 

0.0952 

1.050 

0.8769 

0.0834 

1.038 

0.143 

0.0585 

0.0753 

1.040 

0.8980 

0.0544 

1.032 

0.112 

0.C5CS 

0.0642 

1.030 

0.9281 

0.0420 

1.024 

0.0830 

0.0420 

0.0511 

1.020 

0.9501 

0.0298 

1.0172 

C.0535 

0.0309 

0.0356 

1.010 

0.9744 

0.0139 

1.0065 

0.0255 

0.0175 

0.0198 

1.005 

0.9856 

0.00B7 

1.0042 

0.0122 

0.0025 

0.0097 

Equations  2.48  and  2.49  would  yield  ))  as  a  function  of  (p,  -  Pi)/p4,  if  the  shape  factor  f 
were  known  as  a  function  of  n.  In  the  high-pressure  region,  however,  we  have  already  de¬ 
termined  fl  as  described  above.  Here  Eqs.  2.48  and  2.49  can  be  used  to  calculate  the  shape 
factor.  Ttie  result  is  shown  in  Fig.  2.4.  Also,  for  the  low-pressure  range,  the  asymptotic 
relations  in  Eq.  2.44  can  be  used  to  find  the  shape  factor.  The  value  of  I  in  the  intermediate 
range  can  then  be  approximated  by  interpolation. 

In  the  actual  calculations,  we  used  a  combination  of  the  two  methods  described  above. 
First,  the  integration  constant  c  in  the  asymptotic  expression,  Eq.  2.44,  was  approximately 
determined  by  merging  the  two  parts  of  the  0  curve  graphically,  as  shown  In  Fig.  2.3.  Once 
this  integration  cemitant  is  known,  the  shape  factor  can  be  calculated  for  the  low  pressures 
and  Interpolated  in  the  intermediate  range.  Now  Eqs.  2.48  and  2.49  can  be  used  to  find  0  and, 
subsequently,  an  Improved  integration  constant  c.  The  procedure  is  repeated  until  satisfactory 
convergence  is  obtained.  It  turns  out  that  the  shape  factor  has  a  maximum  and  a  minimum  In 
the  range  of  our  interest  and  varies  between  1.65  and  0.53  (see  Fig.  2.4).  These  calculations 
yield  0  over  the  entire  range  of  our  project.  Subsequently  q|  can  be  determined  by  integration. 
This  concludes  the  major  part  of  the  shock-wave  calculations  since  now  the  peak  pressure- 
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Fig.  2.4 — Stupe  factor,  icduced  toul  energy,  and  time  factor. 


distance  curve  as  well  as  other  shock  parameters  can  be  readily  computed,  as  will  be  de¬ 
scribed  In  the  next  section. 


2.11  SUMMARY  AND  RESULTS 


For  the  benefit  of  the  reader  who  is  less  Interested  in  mathematical  developments,  the 
salient  points  of  the  analysis  are  summarized  here. 

The  main  objective  Is  to  calculate  the  shock-wave  peak  pi  essure  as  a  function  of  distance 
for  ■'  '^.olnt  explosion.  A  convenient  expression  for  the  pressure-distance  relation  may  be  ob¬ 
tained  from  the  fact  that  the  total  energy  of  the  spherical  disturbance  caused  by  the  explosion 
must  be  constant  and  equal  to  Q,  the  hydrodynamic  yield  of  the  explosion.  The  mathematical 
expression  of  this  statement  is 

Q  =  I  ’tJp,  Vi  =  constant,  (2.51) 

where  Pi  =  shock-wave  peak  pressure 
Pi  -  peak  density 

pg  -  density  of  the  undisturbed  water 

r|  -  distance  of  the  shock  front  from  the  center  of  the  explosion 
Vi  =  reduced  total  energy  as  defined  below  (see  Eq.  2.53) 

According  to  the  Ranklne-Hugoniot  conditions,  the  shock  pressure  P|  Is  a  function  of 
p|/pii  and  Is  given  by  the  Ranklne-Hugoniot  adiabatic,  Eq.  1.1.  In  the  thermodynamic  part  of 
the  calculations  this  function  has  been  numerically  determined  and  is  presented  In  Table  3,1. 
ft  will  be  noted  that  P|  varies  between  4po  and  pg.  The  latter  value  correspondc  to  pj 0;  the 
first,  to  P|  =  «;  hence  to  the  instant  of  explosion:  t  =  0  and  ri  =  0. 

It  was  found  advantageous  In  the  analysis  to  use  the  magnitude  p|/pg  as  the  principal  In¬ 
dependent  variable.  (For  simplicity,  we  write  pi  instead  of  the  above  ratio,  assuming  pg  Is  a 
constant  and  known  magnitude.)  A  glance  at  Table  3.1  will  show  that  It  Is  possible,  once  pj  Is 
specified,  to  determine  any  other  variable  of  Interest,  such  as  pressure,  radius,  or  time. 

Tne  reduced  total  energy  ii  depends  on  the  shock  strength;  hence 

Vi  ~  Vi(Pt)-  (2.52) 


Once  this  relation  between  iji  and  p,  Is  found,  the  corresponding  distance  of  tr.e  shock  frent  rg 
Is  easily  found  from  the  energy  equation,  Eq.  2.51, 
t)i  Is  defined  by 


7)1  =  3 


P.-P«Pij\r(/  \r,/ 


(2.53) 


^  3 


\  Pl~  Po 


where  U|  Is  the  particle  velocity  directly  behind  the  shock  front;  the  symbols  without  sub¬ 
scripts  refer  to  points  between  the  front  and  the  center;  <p  Is  the  reduced  veiuulty,  x  Is  the 
reduced  density;  4  Is  the  reduced  pressure;  and  (  Is  the  reduced  distance.  These  reduced 
variables  take  values  only  between  zero  and  unity.  Figure  2.1  shows  the  behavior  of  these 
magnitudes  for  various  values  of  P|.  J  Is  the  reduced  Internal  energy  discussed  In  Sec.  1.2. 

It  Is  a  function  of  three  variables,  namely,  Pi,  it,  and  x<  These  data  are  directly  obtained  by 
the  thermodynamic  calculations  outlined  In  Chap.  1  and  are  shown  in  Fig.  1.4.  Obviously  the 
material  ptesented  in  Figs,  i.4  and  2.1  provides  all  the  Information  necessary  to  calculate  Vt 
by  means  of  Eq.  2.53.  However,  to  obtain  qi,  x,  and  ^  as  functions  of  {  and  pi,  one  must  solve 
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the  three  partial  ditferentlal  equations  of  the  fluid  motion.  In  NAVORO  Report  4164  several 
methods  are  proposed  which  provide  approximate  solutions  of  reasonable  accuracy. 

One  of  the  Interesting  features  is  that  simple  equations  exist  for  the  inclination  of  the  ip, 
and  il  curves  at  {  =  1  (Eqs.  :!.16  to  2.18).  This  means  that  the  behavior  of  the  reduced 
variables  near  {  =  1  can  be  found  In  a  relatively  simple  manner,  which  Is  most  Important  for 
the  finding  of  r/f  by  means  of  the  Integral,  Eq.  2,53,  Since  the  lector  occurs  in  the  integrand, 
values  tor  small  (  contribute  little  to  the  total  value  of  the  integral  and  hence  do  not  need  to  be 
known  precisely.  (This  holds  for  high  pressures  only.  In  the  Intermediate-  and  low-pressure 
ranges  J  increases  lapldly  with  decreasing  Hence,  other  methods  are  used  in  these  ranges.) 
For  greater  accuracy  in  the  evaluation  of  i)|  in  the  high-pressure  range,  we  have  derived  series 
for  the  reduced  variables  which  describe  their  behavior  around  {  =  0.  Combination  of  both 
treatments  by  means  of  polynomials,  which  incorporate  the  information  available  at  {  -  0  as  well 
as  at  t  =  1,  resulted  in  the  expressions  In  Eq.  2.27.  Their  accuracy  is  excellent  in  the  one 
Instance  where  a  comparison  with  the  rigorous  solution  was  possible. 

However,  there  is  one  difficulty  <1  these  otherwise  simple  calculations;  they  involve*  an 
unknown  parameter,  namely, 

^  ,  {--i  -  P»  dill  =  dlniii  ,2  5^, 

d£l 

P)  Pt 

This  magnitude  is  the  inclination  of  the  curve,  when  plotted  vs  (f,-  Po)/Po  in  n  logarithmic 
scale.  Figure  2.2  shows  such  a  plot. 

Actually  ft  is  the  magnitude  of  principal  concern,  because,  once  is  known,  the  im¬ 
portant  magnitude  >)|  can  be  readily  obtained  by  integration.  Therefore  *he  greatest  part  of  the 
analysis  Is  devoted  to  the  task  of  finding  /).  A  different  method  was  used  in  each  of  the  three 
pressure  regions  considered. 

The  hign-prersure  region  corresponds  to  pj  between  4pt  and  1.4pg  and  to  pressures  between 
«  and  484,000  pel.  Here,  nt  «ss  calculated  from  Iq.  3.53  with  various  arbitrarily  assumed 
values  of  0.  Then  Vi  H  piotted  vs  (p,  -  pnVpt  and  an  tj|  curve  in  constructed  which  has  an  incli¬ 
nation  equal  to  the  assumed  0.  Along  this  curve  the  condition  in  Eq.  2.54  is  satisfied,  and  this 
0  is  the  solution  of  our  problsm.  The  interpolation  iiieihod  used  to  select  0  Is  Illustrated  in 
Fig.  2.3. 

The  low-preseure  region  extends  from  pi  l.Olpg  down  to  pi  =  p«  and  corresponds  to 
pressures  between  3260  psl  and  zero.  For  this  region,  the  following  asymptotic  relation  can  be 
derived: 


This  determines  0  .ts  a  function  of  p^  and  an  integration  constant  c  which  is  to  be  determined, 
figure  2.3  shows  the  form  Jt  the  0  curve  when  plotted  vs  the  logarithm  of  (pj  -  Po)/pO'  A  change 
ill  c  moves  this  curve  horizontally  either  to  the  left  or  to  the  right  without  changing  the  ordl- 
na'-is.  This  shows  how  the  0  curve  can  be  approximately  drawn  for  both  the  low-presrure  and 
the  lutermediaie-piessuie  ranges;  The  asymptotic  pressure  curve  ts  moved  into  such  a  posl- 
tic-*  that  a  emooth  curve  is  (Stained  from  the  high-pressure  down  to  the  low-pressure  range. 
The  curve  thv.s  obtained  corresponds  to  the  fine  line  in  Fig.  2.3.  (TTie  figure  actually  shows  not 
this  crude  Intermediate  approximation  but  the  final  result  obtained  with  the  use  of  the  subse¬ 
quently  described  method.) 

The  intermediate-pressure  range  has  been  calculated  by  meane  of  a  "peak  approximation." 
With  this  term,  we  refer  to  the  theories  in  which  an  attempt  is  made  to  describe  the  shock- 
wave  propagation  Ir  termr.  of  shock  parameters,  such  as  peak  pressure,  time  constant,  or  re- 


*Compare  Eqs.  2.16  to  2.18,  2.10,  and  2.14  of  this  report. 
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iaied  magnitudes.  This  is  od'y  possible  I!  assumptions  as  to  the  shape  of  the  wave  profile  are 
made,  and  they  are  aopruximatloas  fur  that  reason.  Examples  are  the  highly  successful 
Kirkwood -Bcthe  and  Kirkwood- Brinkley  theories.  We  have  used  In  this  project  the  more  recent 
theory  of  Snay  and  Matthias.'  Adapted  to  oi<r  problem,  this  theory  yields: 


t>  =  - 


W,  +  W2OIX 


Pi  -  P*  4 


d  In 
d  In 


n  ax  - 1  /  w.  \ 

Pi-Po  Wj+WjffxV  *  orx  •  y 


Co 


(2.S8) 


where  Wj  to  Wj  -  functions  of  pj,  listed  la  Table  2.1 
ox  =  time  factor  -  Ff/c^S 
Cj  =  sound  velocity  at  P|  =  0 
e  =  time  constant  =  p  -  (8p/6t) 

f  -  shape  f.ictor,  see  Eq.  2.45  (f  Is  unity  for  an  exponential  wave,  xero  for  a 
triangular  wave) 

riiK  -  fd  In  P|)/<d  In  p,),  listed  in  Table  3.1. 

At  the  end  of  the  high -pressure  range,  Pi/po  -  1.8,  we  know  all  three  magnitudes  involved  In 
Eqs.  2. 56,  namely  fi,  ax,  and  f.  The  latter  two  are  the  Initial  conditions  for  the  above  differen¬ 
tial  equation.  To  Integrate  this  equation,  f  must  be  known  as  a  function  of  p,.  Assumption  of  a 
constant  f  would  have  been  a  good  approximation.  We  went  even  further  by  estimating  f  as  a 
function  of  p,,  as  described  In  Sec.  2.10  and  shown  in  Fig.  2.4.  Finally,  integration  of  Eqs. 

2.56  yields  0  as  well  as  ax  from  which  the  time  constant  B  can  be  obtained.  The  results  are 
given  in  'fable  3.1. 

Subsequent  straightforward  calculations  yield  Pi,  Eq.  2.54,  and  the  shock-front  distance  f), 
Eq.  2.51.  This  yields  the  desired  pressure-distance  relation.  The  other  shock-wave  param¬ 
eters,  such  as  time  constant  and  shock-wave  energy,  can  now  be  calculated  without  difficulty. 
This  will  be  explained  In  the  following  paragraphs,  together  with  the  presentation  ot  the 
numerical  results. 

The  peak  pressure -distance  curve  for  a  point  explosion  with  a  yield  of  30  kt  of  TNT  Is 
presented  l'><  Fig.  2.5  together  with  the  experimental  evidence  obtained  In  Operation  Wigwam, 
as  given  In  the  preliminary  version  of  reference  15.  Also  shown  are  the  results  of  the  pre¬ 
liminary  calculations  submitted  prior  to  the  test  as  "predictions.”  The  two  curves  differ  In  the 
range  of  Interest  by  very  minor  amounts,  and  either  one  is  in  good  agreement  with  the  experi¬ 
mental  results. 

The  pressure-distance  relation  for  any  other  yield  Q  can  be  easily  calculated  from  the 
values  of  r|/Q'’  which,  together  with  the  peak  pressure,  are  given  In  Table  3.1.  For  low  pres¬ 
sures  (3000  psl  and  below),  this  relation  cxn  be  approximated  as  follows: 


Pi 


=  4.608  X 


--  18,800 


(2.57a) 

(2.57b) 


where  P|  Is  In  pounds  per  square  Inch;  r|  is  in  feet;  Q  is  In  kllotons  of  TNT,  where  1  kt  Is  lO" 
g-cal  (~4.20  \  10”  ergs);  and  W  Is  tn  pounds  of  TNT  (1  kt  =  2.205  x  10*  lb).»  Comparison  of 

*A  detonation  energy  ot  lu’  cal  g  hes  been  assumed  for  TNT.  Previously  published  values 
which  are  up  to  10  per  cent  higher  than  this  are  based  on  theoretical  calculations  which  are 
superceded  todav.  A  reliable  value  is  not  known,  and  there  are  indications  that  even  10*  cal/g 
Is  too  high.  'This  shows  cleat  ly  that  the  choice  of  kllotons  of  TNT  as  an  energy  unit  Is  :m  un¬ 
fortunate  one. 
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Eq.  2,57b  with  the  cur  responding  formula  for  the  shock-wave  peak  pressure  produced  by  a  TNT 
explosiun, 


p,  =  21,600 


(2.58) 


shows  that  for  a  point  explosion  an  energy  equivalent  to  that  of  1.446  lb  of  TNT  la  necessary 
to  produce  the  same  shock  pressure  as  does  1  lb  of  exploding  TNT.  Or,  in  other  words,  the 
shock-wave  peak  pressure  of  a  point  explosion  with  a  yield  of  1  kt  of  TNT  Is  equal  to  that  of 
an  exploding  0.692-kt  TNT  charge. 

The  shock-wave  peak  pressures  observed  at  Wigwam  correspond  to  those  produced  by  a 
TNT  explosion  of  46.2  x  10*  lb  charge  weight.**  The  actual  yield  of  the  Wigwam  explosion  is 
not  accurately  known.  The  followiiig  values  were  quoted  at  the  Wigwam  Project  Officers’ 
Conference  at  the  Naval  Radiological  Defense  Laboratory,  San  Francisco,  in  October  1055: 

Radiochemical  yield: 

Los  Alamos  Scientific  Laboratoi y,  32  x  3,2  ki  of  TNT 

Naval  Research  Laboratory,  35  -  2.5  +  10  kt  of  TNT 

Hydrodynamic  ylela: 

Armour  Research  Foundation,  30.5  x  1  kt  of  TNT. 

Using  the  LASL  value,  it  is  found  that  according  to  the  experimental  evidence  of  Operation 
Wigwam,  an  atomic  underwater  explosion  with  a  yield  of  1  kt  of  TNT  would  produce  the  same 
shock-wave  peak  pressure  as  a  TNT  explosion  of  0.65  ±  0.13  kt  charge  wel^t.  The  calculated 
value  (0.692  kt)  is  within  the  limits  of  the  experimental  error.  This  error  appears  to  be 
large,  because  expressing  the  shock-wave  peak  pressure  in  terms  of  energy  amplifies  the 
errors  of  the  pressure  measurements  by  the  power  3/1.13  =  3.555.  For  the  same  reason  it  Is 
more  difficult  to  predict  accurately  the  energy  equivalent  than  the  shock-wave  pressure. 

If  the  calculated  curves  in  Fig.  2.5  were  drawn  for  a  yield  of  33  kt  of  TNT  Instead  of  30  kt 
of  TNT.  the  pressures  would  b-;  2.5  per  cent  higher — a  change  which  would  hardly  affect  the 
agreement  with  the  experimental  points  shown.  Since  the  question  of  the  actual  yield  was  not 
settled  at  the  time  of  this  writing,  all  data  In  this  report  are  presented  for  30  kt  of  TNT,  the 
yield  for  which  the  calculations  were  originally  made. 

The  time  constant  Is  defined  by 


Calculated  numerical  values  i-f  0/Q'’  as  well  as  $  for  30  kt  of  TNT  are  listed  In  Tabic  3.1. 
Figure  2.6  shows  a  c.':mparlson  with  the  test  results:  the  calculated  time  constant  is  about  10 
per  cent  too  high.  Since  most  shock-wave  theories  give  poor  results  for  the  time  constant,  the 
agreement  may  be  co.asldered  satisfaciory. 

The  arrival  lime  for  the  shuck  front  at  a  point  at  a  distance  r|  is  given  by 


t  r 
•I  u 


(2.60) 


where  U  Is  the  propagation  velocity  of  the  shock  front 


U  = 


(Pi  Pi  V* 
VPo  Pi  -  PiJ 


(2.81) 


Introduction  of  the  variables  11  and  m  yields 


Q’’  1 

3(po)'>  (V)’’ 


i  (^)'(>h)’>(P.)’‘''”^ 


dp,. 
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(2.62) 


o  o  o  o  o  o 

O  O  «  lO  CVJ  — 

oasn)  iNvxsNoo  bnu 
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Numerical  rcsulls  for  l  as  well  as  t  q’’  arc  slven  In  Table  3.1.  For  high  pressures  and  small 
distances,  the  arrival  turn  i,..ii  be  obtained  by  means  of  the  Taylor  solution  (reference  13,  Eqs. 
7  and  8).  In  our  notation. 


t  ■ 


2/5 

lA 


=  £.588  X  10"* 


=  0.830 


(2.63) 


where  the  numerical  values  apply  to  p|  =  4pj.  t  in  milliseconds,  Q  in  kilotons  of  TNT,  W  In 
pounds  of  TNT,  and  r,  In  feet. 

For  moderate  and  large  distances  Eq.  2.62  gives  good  agreement  with  the  test  results  re¬ 
ported  In  the  preliminary  version  of  reference  16  (see  Fig.  2.7).  However,  for  close-in  dis¬ 
tances  both  Eqs.  2.62  and  2.63  '»rc  In  poor  agreement  with  the  test  results  (Fig,  2,7).  This 
disagreement  Is  probably  due  to  the  Idealizing  assumption  of  a  point  explosion  with  an  Isentroplc 
exponent  of  Vj.  For  real  explosions  this  value  may  be  higher. 

The  shock-wave  energy  flux  (commonly  called  shock-wave  energy)  is,  for  low  pressures, 


esw  =  —  fp^dt. 


(2.64) 


This  magnitude  Is  equivalent  to  the  Integral  Qg  In  Eq.  2.37  and  can  be  expressed  as  follows: 


-  .  VS  Q'‘  (Q'>\ 

•  7,  4T  \1-) 


(2.63) 


where  qsls  the  reduceti  shock-wave  energy  defined  by  Eq.  2.38.  The  ratio  (l/4»)(7)s/»)|)  Is 
listed  In  Table  3.1  It  Is  a  slowly  varying  function  of  p|  or,  what  amounts  to  the  same  thing,  of 
Q'l.'r,.  For  pressures  below  3000  psi  we  have  approximately 


—  ^  =  2.15  X  )  . 

4x  ti,  \  ri  / 

and,  by  combination  with  Eq.  2.65, 

E^i\^•  =  5.522  X  10*  q'> 

=  1963W''  (^)  , 


(2.66) 

(2.67a) 

(2.67b) 


wher"  *'eiV  is  in  lnch-)ioundB  pe;  square  Inch,  rj  Is  In  feet,  Q  is  In  kilotons  of  TNT,  and  W  Is  In 
pounds  of  TNT.  It  should  be  noted  that  these  theoretical  .’orraulas  apply  to  the  total  shock-wave 
energy.  In  the  empirical  formula  for  TNT, 


'SW 


=•  2410W 


(2.68) 


the  Integration  In  Eq.  2.64  Is  carried  up  to  68  and  not  to  Inllnlty.  However,  for  all  pr.nctlcal 
purposes  this  is  equivalent  to  the  total  shock-wave  energy.  The  experimental  values  shown  In 
Fig.  2.8  are  obtained  by  Integrating  up  to  the  time  of  the  surface  cutoff  only;  this  occurred  at 
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alK)ut  I  ■■  ie  io  \  -  30.  They  represent  only  fractions  of  the  total  energy,  although  very  large 
fractions.  Considering  these  factors  of  uncertainty,  the  agreement  between  theory  and  experl- 
n.ent  Is  fairly  good. 

An  explosion  of  4S.76  x  10*  lb  of  TNT,  which  would  give  the  same  shock-wave  peak  pressure 
as  a  30-kt  TNT  point  explosion,  would  not  produce  the  same  shock-wave  energy  (see  Fig.  3.8). 

An  explosion  of  54.06  x  10*  lb  of  TNT  would  be  necessary  to  Jo  this,  as  can  be  found  from  £qs. 
2.68  and  2.67h. 
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CHAPTER  3 


BUBBLE  PHENOMENA 


3.1  INTRODUCTION 

The  nature  of  the  bubble  formation  by  a  point  explosion  differs  in  some  respects  from  that 
of  conventional  explosives.  In  the  latter  case  the  reaction  products  of  the  explosion  expand  and 
push  the  water  aside,  thus  forming  a  gas  sphere  which  then  performs  the  well-known  bubble 
pulsations.  The  bubble  Interface  always  consists  of  the  same  particles,  namely,  those  whicii 
were  In  contact  with  Ihe  charge  before  the  explosion.  Vaporisation  of  the  water  Is  entirely 
unip'p''rtant. 

In  the  case  of  a  point  explosion,  the  water  Is  pushed  aside  by  expanding  steam.  The  heat 
which  vaporizes  the  water  stems  from  the  energy  dlsnipated  at  the  shock  front  and  Is  that 
energy  which  the  shock  wave  leaves  behind  In  the  form  of  thermal  energy  after  the  medium  has 
reexpanded  to  the  Initial  pressure.  Consequently,  the  water  has  a  higher  temperature  In  this 
state  than  it  had  before  the  passage  of  the  siiock.  This  temperature  Increment  decreases  with 
Increasing  dlslance  from  the  center.  Near  the  center  It  is  high  enough  not  only  to  vaporize  but 
actually  to  decompose  the  water.  At  greater  distances  this  temperature  increment  eventually 
drops  to  a  value  which  corresponds  to  the  boiling  point  of  water.  This  condition  determines  the 
bubble  radius.  Since  the  boiling  point  Is  a  function  of  pressure  as  well  as  temperature,  the 
mass  of  water  which  Is  evaporated  depends  on  the  prevailing  pressure,  which  changes  during 
the  expansion  of  the  bubble.  Thus  the  Interface  is  not  formed  by  the  same  layer  of  water  but  is 
transferred  from  one  set  of  particles  to  another.  At  the  moment  of  the  bubble  maximum  the 
greatest  mass  of  water  is  In  vapor  form.  The  radius  of  this  maximum  bubble  Is  the  parameter 
of  primary  interest.  All  the  other  bubble  parameters,  such  as  the  first  period  and  the  bubble 
energy,  as  well  as  the  numerical  values  of  the  energy  partition,  are  readily  deduced  once  this 
magnitude  Is  found. 


3.2  DETERfvQNATICW  OF  MAXIMUM  RADIUS 

The  energy  which  has  been  Imparted  to  the  water  by  the  explosion  can  he  subdivided  into 
three  portions; 

1,  Shock-wave  energy,  Qsw. 

2,  Dissipated  energy,  Qduf. 

3,  Bubble  energy,  Og, 

The  bubble  energy  Is  most  coitveitlemly  found  from  the  dlflercnce  of  total  energy  and  the 
first  two  energy  terms  noted  above: 

Qii  “  Q  -  Qsw  -  Qdifi-  (S«l) 


sa 
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The  shock-wave  energy  fraction  Is,  according  to  Eqs.  2.38  and  2.3, 


5siSi  a  ^ 
<5  JJi 


(3.2) 


This  sr.sgnltu.'ic  is  a  function  of  distance  and  has  been  calculated  in  the  previous  chapter. 
Numerical  results  can  be  obtained  from  Table  3.1. 

According  to  Eq.  2.37,  the  dissipated  energy  for  a  shock  winch  has  traveled  from  the 
distance  r*  to  r*  is 


=  <jp,  jTj*  hr' dr;, 


(3.3) 


where  h  Is  the  dissipated  enthalpy  Increment,  Eq.  1.14,  and  r|  is  the  distance  of  the  shock  front 
from  the  center  of  the  explosion.  With  the  use  of  Eqs.  1.12,  2.2,  and  2.34,  we  obtain  for  the 
dissipated  energy  fraction 


(3.4) 


The  lower  limits  of  the  integrals  in  Eqs.  3.3  and  3.4  refer  to  the  point  at  which  the  shocjc-wave 
energy,  Eq.  3.1,  is  calculated.  The  upper  limits  of  these  integrals  must  be  chosen  Ui  such  a 
way  that  the  Integral  covers  only  those  particles  which  are  in  the  liquid  state  at  the  moment  of 
the  bubble  maximum,  since  any  particle  which  vaporises  belongs  to  the  interior  of  the  bubble 
and  its  energy  is  counted  as  bubble  energy.  Therefore  r*  designates  the  layer  of  particles 
which  forms  the  in'eiface  of  the  bubble  at  Its  maximum  expansion,  and  the  sphere  harin';  this 
radius  contains  the  mass  of  water  evaporated.  When  the  shock  front  has  traveled  to  the 
distance  r],  it  encompasses  a  sphere  of  average  density  p,,  the  normal  density  of  the  water 
before  the  explosion.  Later,  this  sphere,  consisting  of  the  same  particles,  has  expanded  from 
the  radius  r*  to  the  maximum  bubble  radius  Am.  Of  course,  the  mass  of  these  two  spheres  is 
the  s,ame,  and  therefore 


m  ”  dm  ■  4x^5  ^  pr'dr 

r*  « 

'**  V  P**"'  ‘‘‘■i  • 


The  volume  of  the  bubble  at  its  maximum  is  then 


<»  /o'  Pov(r)r]dr,. 


(3.8) 


where  v  is  the  specific  volume  of  the  steam  Inside  the  bubble  at  the  moment  of  the  maximum 
expansion.  We  assume  that  at  this  moment  the  pressure  P],(  is  constant  throughout  the  bubble. 
Since  the  bubble  interface  consists  of  saturated  liquid,  this  pressure  is  the  saturation  pressure 
of  water  in  the  thermodynamic  state  prevailing  at  the  interface.  The  pressure  Pg  and  the 
specific  volume  vs  of  the  saturated  liquid  are  functions  of  the  entropy  S  only  and  can  be  readily 
obtained  from  the  steam  tables,  such  as  In  reference  17.  Thus 


Pm  =  Ps(Sa„) 
v(Am)  =  vs(Sa„). 


(3.7) 
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Also  liste'l  ill  tho  <;ir>am  tables  is  the  c>peci{tc  volume  ol  steam  at  constant  pressure  as  a  func¬ 
tion  of  entropy.  Therefore  we  have 


v(r)  =  vlS(r),P;^J 


(3.8) 


---  vlS(r),S(AM)I. 

Since  the  entropy  is  constant  alosg  the  patti  ct  a  particle,  once  the  snocK  has  passed  over  It,  we 
have 


S(Am)  =  C{r\) 

S{r)  =  S(r,)  (3.9) 


and 


v(r)  =  vlS(r,),S(r;)|. 

This  shows  that  v(r)  can  be  determined  with  the  use  ol  the  steam  tables  and  Table  11.1.  (The 
steam  tablss  give  tne  entropy  increment  above  0°C,  whereas  Table  3.1  gives  that  above  8°C. 
Therefore  S(8°C)  -  S(0°C)  must  be  added  to  the  latter  values  to  bring  them  on  an  equal  basis.) 

The  exact  evaluation  of  the  integral  In  Eq.  3.6  is  difficult  because  v  increases  to  infinity 
when  rj  approaches  zero.  A  rather  crude  approYtmution  was  employed,  using  Simpson’s  rule 
for  the  evaluation  of  this  integral.  The  result  ts  shown  in  Fig.  3.2  (the  curve  marked  Eq.  3.6), 
where  the  dimensionless  magnitude  A^/r^  is  plotted  vs  p/po. 

To  obtain  a  solution  of  our  problem,  we  need  another  independent  expression  foi  the 
magnitude  A,vi/r*.  According  to  basic  bubble  theory  (Eq.  8.3  of  reference  14),  the  maximum 
bubble  radius  and  ‘he  bubble  energy  are  related  by 


A*  -  —  rt 


(3.10) 


where  Po  is  the  absolute  hydrostatic  pressure  at  firing  depth  and  is  a  fafctor  which  accounts 
for  the  internal  energy  within  tne  bubble.  (In  the  bubble  theory,  absolute  energies  arc  usually 
used,  whereas  the  cncrgle.'>  Q,  Qsw>  imd  Qsare  excess  energies.  These  must  also  be  accounted 
for  by  a  .) 

Introductior.  of  Eq.  2.2  into  Eq.  3.10  yields 


=  16,284a 


>7iPi  Pi  -  Pd 

M  -fT-TT- 


(3.11) 


where  Zg  Is  the  total  hydrostatic  head  In  feet  and  P|  is  in  kilobars.  Figure  3.2  (the  curve 
marked  Eq.  3.11)  shows  the  resuito  obtained  from  this  equation  together  w'th  Eqs.  3.1  to  3.3. 
The  curve  holds  for  Zg  -  ?,033  ft  and  wjj  -  0.8.  The  lower  limit  of  the  integral  in  3q.  3.4  was 
oet  at  p?  •••  l.Olpg.  From  Table  3.1,  one  finds  the  corresponding  shock-wave  energy  fraction 
Qsw 'Q  -  0-20. 

The  two  curves  for  A'.t/i  f  Intersect  at  tne  point 


£i  =•  1.597  and  =  6.87.  (3.12) 

Po  1 

This  establishes  the  solution  of  our  problem  and  dnormines  immediately  the  following  magni¬ 
tudes,  as  illustrated  in  Figs.  3.1  and  3.2: 
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Fig.  3.? — Detecinlnallon  af  the  maximum  bubble  ndlui. 
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(3.13) 


S!Am)  =  0.415  +  S  (O'C) 

^  =  0.390. 

w 

3.3  BUBBLE  PARAMETERS  AND  ENERGY  ?ARTJ'nON 

From  the  results  !e  Fqs.  3.12  and  3.13  Wf  Itiid  the  maximum  bubble  radius  lor  Q  "  30  ht  of 
TVT  und  7,  =  2033  .’I: 


Am  =  373  ft  ;S.H) 

fhe  corrcs.xinding  period  of  the  first  bubble  pulsation  T  is  readily  obtained  from  the  ex¬ 
perimental  evidence  that  the  rntlo  of  the  period  constant  K  and  the  radius  constant  J  are  al¬ 
most  the  same  (or  all  high  explosives  ter.ted  so  fai.'*  Using  the  value  C.34S,  which  holds  for 
TNT,  we  have 


K  =  C.345J 


T(7j^_ 

W'j 


0.345 


Am(Z|.)'- 

WC~ 


(3.15) 


T  =  0.345 


Am 

(Z,)'?' 


Tils  yields  the  first  bubble  period  for  a  30-kl  TNT  point  exoloslon  at  a  depth  of  2000  ft? 


T  =  2.89  sec. 


(3.1 5a) 


Here  the  question  arises  whetner  corrections  to  the  period  lor  r’jrface  or  bottom  effects 
should  be  made.  (The  maximum  radius  is  not  changed  by  these  effects.)  An  approximate 
formula  for  sue!'  a  correction  la  (see  Eq.  18,  p.  5-3  of  reference  19): 

T  -  T,  ^1  -  0.2  ^  +  0,2  .  (3.16) 

wiiere  T,  -  free  water  period  as  calculated  above 
=  maximum  bubble  radius 
D  =  depth  of  explosion  below  water  surface 
B  :>  depth  of  bottom  below  center  of  explosion. 

This  and  similar  formulas  found  In  the  literaturu  btuiry  overcorrect  the  bottom  effect. 

They  also  overcorrect  the  surface  effect  for  large  high-exploeivo  charges  once  the  bubble  It 
several  maximum  radii  below  the  surface.  The  latter  rlfect  >s  not  eo  well  established  sc  the 
first  one  and  Is  not  mentioned,  in  the  literature,  but  It  is  evident  from  Fi|;.  8.21  of  refe'ence  14. 
Fcr  configurations  similar  to  that  considered  here  (depth  about  sVi  maximum  radii),  the  un- 
corrected  equation  seems  to  give  tlie  more  accurate  result.  The  period  corrected  for  bottom 
and  surface  effects  would  be  2.79  sec,  i.e.,  about  1  per  cent  smeller. 

In  Sec.  3.2  the  partition  of  the  three  energy  terms  has  been  calculated;  see  also  fig.  3.1. 
We  summarise: 
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(•'pergy  balance: 

Bubble  enerK>'.  39  per  cent  of  tin  total  energy 
Dissipated  energy,  41  per  cent  ot  the  total  energy 
Shock-wave  energy,  20  per  cent  ot  the  total  energy. 

The  last  two  terms  deper.d  on  the  distance  to  which  the  shock  has  traveled.  They  hold  for 
b|  =  l.Olpo  (which  is  the  value  of  the  lower  limit  used  in  the  calculation  of  the  Integral,  Eq.  3.4). 
This  corresponds  to  the  shock-front  distance 

r,  =  812.4Q'> 


or 


r,  =  1903  ft  for  n  30-kt  TNT  point  explosion. 

At  larger  distances  the  shock-wave  energy  decre.'tses  as  shown  in  Table  3.1,  whereas  the 
dlseipatcd  energy  increases  correspondingly. 


3.4  BUBBf  .E  PRESSURE  AND  TEMPERATURE,  MASS  OF  WATER  VAPORIZED 

In  Sec.  3.2,  we  have  found  the  entropy  on  the  bubble  interface  at  the  moment  of  the  maxi¬ 
mum  expansion.  From  the  steam  tobies  we  can  find  the  corresponding  pressure  and  tempera¬ 
ture  of  the  saturated  liquid.  The  llrst  magnitude  is  equal  to  the  bubble  pressure  at  the  moment 
of  the  bubble  maximum,  and  the  latter  r  •era  to  the  temperature  on  the  bubble  interface  at  ilic 
same  moment.  From  the  jteam  tables  we  find 


Pm  •-  52.4  psia 
T(Am)  =  7.84'F. 


(3.17) 


At  the  maximum  bubble  radius  the  bubble  pressure  has  therefore  dropped  to  '/|;j  of  the 
hydrostatic  pressure  and  is  constant  throughout  the  inside  of  the  b>  b'.>lc.  The  temperature  is 
constant  only  within  the  shell,  adjacent  to  the  bubble  Interface,  which  contains  moist  steam. 
From  this  point  on,  tho  temperature  increases  rapidly  with  decreasing  distance  from  the 
center  and,  theoretically,  reaches  Infinitely  high  values  at  the  center. 

The  density  of  the  medium  inside  the  oubbie  is,  at  the  Interface,  that  of  the  saturated  liquid 
(about  0.925  g/cc).  The  density  decreases  rapidly  with  deceasing  distance  and  vanlshiss  at  the 
centir. 

Wiien  the  shock  fruni  has  traveled  to  the  distance  rf,  it  has  reached  the  layer  of  particles 
which  will  be  on  the  bubble  interface  at  the  moment  of  the  maximum  expansion.  The  interior 
of  the  sphere  of  radius  r*  r/Ul  be,  therefore,  the  interior  of  the  bubble.  Kence,  the  mass  of 
water  evaporated  is 


m  ^  ^  P,(rr)’- 

We  have  found  in  Sec.  3.2  that,  for  a  firing  depth  of  2000  ft, 
rj  =  17.8Q'’. 

This  corresponds  to  a  radius 

rJ  =  54.7  ft  for  Q  =  30  k(  of  TNT. 

With  a  density  fur  sea  water  of  1.025  g/cc,  we  find 
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Maximum  mass  o(  water  evaporated*  =  19.9  metric  kt. 

It  Is  Interesting  that  the  ball  ot  water  which  Is  evaporated  has  Just  the  slse  of  a  SO-kt  TNI' 
sphere.  (Tne  radius  of  the  latter  Is  5S  ft  for  a  loading  density  of  1.52  g/cc.) 

At  the  moment  of  the  bubble  maximum  the  greatest  amount  of  water  is  evaporated.  When 
the  bubble  contracts,  the  pressure  Increases  and  steam  condenses  at  the  Interface.  Thus  the 
Interface  Is  transferred  to  particles  which  were  previously  inside  the  bubble.  The  condense* 
ilon  will  cease  when  the  bubble  pressure  and  temperature  have  Increased  to  the  critical  point 
of  waior.  Beyond  this  point  no  condensation  is  possible,  and  the  mass  of  the  bubble  remains 
constant  during  any  lurcher  contraction.  This  mass  will  always  remain  in  the  vapor  state  as 
long  as  the  flow  pattern  of  the  bubble  pulsation  Is  Irrotatlonal,  l.e.,  as  long  as  there  is  no 
mixing  between  the  steam  and  the  surrounding  cold  water.  This  mass  of  steam  Is  easily  cal* 
culated  from  the  cntropy-dtstrnce  relation  which  can  be  obtained  from  Table  3.1.  We  must 
determine  that  shock  distance  r|  for  which  the  entropy  behind  the  front  Is  equal  to  the  entropy 
at  the  critical  point.  Ihe  latter  is  Ser  =  1.058  S  (0*C).  The  corresponding  shock  radius  Is 
found  to  be  38  ft,  and  the  mass  of  water  which  theoretically  always  remains  in  the  v^ior  phase 
is  5.87  metric  kt. 

Actually,  there  will  be  a  strong  mixing  of  different  water  a.nd  steam  layers  near  the 
bubble  minimum  for  two  different  reasons.  The  first  Is  the  Instability  of  the  interface  during 
the  period  of  time  when  the  bubble  is  near  its  minimum.  This  produces  the  disintegration  ot 
the  interface  into  a  spray  which  is  thrown  into  the  interior  of  the  bubble.  Obviously,  this  brings 
about  a  thorough  mixing  of  the  steam  with  the  surroundir^f,  cooler  water.  The  second  phenom* 
enon  causing  mixing  is  the  distortion  ot  the  bubble  sluqpe  in  the  gravitational  field.  When  the 
originally  spherical  bubble  contracts,  its  lower  boundary  moves  Inward  faster  than  the  other 
points  ot  Its  surface.  The  cross  section  ot  the  bubble  assumes  the  shape  of  a  kidney,  and  finally 
the  lower  boundary  Impinges  upon  the  upper  boundary,  causing  a  vast  amount  ot  turbulence 
Inside  and  outside  of  the  bubble.  It  Is  for  these  reasons  that  the  concept  of  the  mase  of  water 
which  remains  in  vapor  form  is  misleading  and  that  considerable  condensation  must  be  expected 
at  the  bubble  minima. 

3.5  BUBBLE  MIGRATKXf  AMD  LATER  BUBBLE  PHA8I8 

NAVORO  Report  <185  describes  a  method  for  the  calculation  of  the  bubble  oscUlsUon  and 
migration  for  the  second  and  later  cycles.  The  essential  part  of  these  calculations  is  the  de* 
termination  of  the  bubble  energy  for  the  subsequent  cycle.  At  each  minimum  the  bubble  energy 
is  reduced  owing  to  the  emission  of  the  buUile  pulse  and  an  energy  dissipation  which  is  not 
entirely  understood  at  the  present  time.  Probably  it  is  closely  connected  with  the  phenomena 
described  at  the  end  of  Sec.  3.4,  In  NAVORO  Report  4185  these  energy  losses  are  determined 
by  a  semlemplrlcal  method  which  uses  experimental  data  obtained  with  high  exploslvea.  Obvl* 
ously,  these  calculations  are  not  valid  for  steam  bubbles.  Steam  bubbles  suffer  the  sasfs  energy 
losses  as  gas  bubbles  do,  but  in  addition  thers  is  condensation  which  damps  the  oaclUatlon  even 
more.  Therefore,  the  calculations  tor  gas  bubbles  may  be  considered  as  an  upper  limit  for  the 
periods  and  maximum  radii  of  steam  bubbles. 

The  following  results  were  obtained  for  gas  bubbles: 

Second  cycle: 

Maximum  radius,  368  ft 
Period,  3.40  sec 

Bubble  energy,  84  per  cent  of  the  bubble  energy  of  the  first  cycle 

Migration  between  first  and  second  bubble  maxima,  675  ft  (There  is  no  appreciable 

migration  up  fo  the  moment  ol  the  first  bubble  maximum.) 


•The  bubble  co.itains  a  large  ar.iount  of  "wet”  steam.  Her.ce,  strictly  speaking,  only  a 
fraction  of  this  macs  Is  in  the  vapo.*  phase. 
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Thtrii  rvc'c: 

Maximum  ridius.  330  ft 
Period,  4.38  sec 

Bubble  energy,  18  per  cent  of  the  bubble  energy  of  the  first  cycle 
Migration  between  second  and  third  bubble  maxims,  793  ft. 

The  gas  bubble  reaches  tht  water  surface  shortly  after  the  third  minimum,  about  10.7  sec 
after  the  explosion.  These  results  are  graphically  illustrated  in  Fig.  3.3.  It  will  be  noted  that 
the  periods  of  the  later  cycles  are  increased,  although  the  bubble  energies  are  smaller.  This 
is  because  the  bubble  has  migrated  into  shallow  water  where  the  lower  hydrostatic  pressure 
causes  the  bubble  to  pulsate  more  slowly. 

Very  similar  data  were  .ibtalncd  expertr.'.c.ilai;y  in  the  vacuum  tank.”'**  In  cn?  of  these 
experiments,^'  electric  sparks  were  used  to  deliver  the  energy  of  explosion.  The  bubbles 
produced  in  this  way  are  steam  bubbles  which  should  behave  very  much  like  the  bubbles  from 
atomic  underwater  explosions.  Unfortunately,  condensation  phenomena  are  not  correctly 
scaled  in  such  tests.  To  scale  gravity,  the  pressure  above  the  water  surface  must  be  reduced 
.<>0  far  that  the  vapor  pressure  of  water  is  almost  reached.  However,  for  similarity  of  con¬ 
densation  phenomena,  the  ratio 

Total  pressure  at  firing  depth 
Vapor  pressure 

must  be  the  same  in  full  scale  as  in  the  model  test.  Since  it  is  not  possible  to  reduce  the  vapor 
pressure  of  water  by  tbe  same  amount  as  the  total  hydrostatic  pressure  in  the  tank,  condensa¬ 
tion  cannot  be  scaled.  The  conditions  in  the  vacuum  tank  resemble  explosions  in  almost  boiling 
water,  and  little  condensation  is  expected  under  such  circumstances.  In  fact,  the  test  results 
obtained  with  sparks  and  with  high-explosive  charges  are  similar  and  are  in  good  qualitative 
agreement  with  the  calculated  data  above. 

If  in  a  model  experiment  which  employs  electric  sparks  the  air  pressure  above  the  water 
surface  is  not  reduced,  gravitational  phenomena  are  not  scaled  and  we  have  the  case  of  a  non- 
migrating  buDble.  The  scaling  of  condensation  processes,  however,  is  much  Improved  but  still 
not  perfect.  Such  tests  showed  strong  condensation  effects  in  the  second  and  third  cycles.  The 
later  periods  were  substaiitlally  less  than  those  of  corresponding  high  explosives,  which  indi¬ 
cates  strongly  reduced  bubble  energies.  (The  bubble  energies  of  nonmigrating  bubbles  are 
proportional  to  the  cubes  of  tbe  periods.)  The  pulsations  virtually  ceased  after  the  third  cycle, 
and,  apparently,  most  of  the  vapor  was  condensed. 

The  following  figures  give  a  summary  of  the  results  of  these  tests  as  well  as  data  re¬ 
ferring  to  TNT  explosions: 


Periods  relative  to  the  period  cf  the  first  cycle 


Cycle  1 

Cycle  2 

Cycle  3 

(a)  ’lonmlgratlng  TNT  bubble 

1 

0.72 

0.59 

(b)  Nonmlgratlng  steam  babble 

1 

0.45 

0.21 

(c)  Migrating  TNT  bubble  (Wigwam  conditions) 

1 

1.18 

1.52 

(d)  Crude  estimate  for  migrating  steam 

1 

0.74 

0.54 

bubble:  (d)  =  (c)  x  ((b)/(a)| 


sources:  (a)  reference  22,  (b)  preliminary  evaluation  of  Hudson's  teats  described  in  refer¬ 
ence  21,  (c)  calculations  based  on  NAVORD  Report  4185  as  described  in  this  section. 

The  first  period  observed  in  Operation  Wigwam  was  2.87  sec,  the  second,  2.6  sec,  and  the 
third,  1.9  sec  (reference  15).  The  ratios  of  these  latter  periods  to  the  first  period  are  C.91  and 
0.66,  respectively.  These  values  are  higher  than  the  above-mentioned  cnide  estimates  fur  the 
migrating  steam  bubble  but  considerably  lower  than  those  for  the  migrating  gas  bubble.  Tnis 
indi.'ttes  that  substantial  rondcn.salion  must  have  occurred  in  Wigwam,  altbnugh  not  quite  so 
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!.  ii‘'h  as  In  tlic  case  of  Hudssn’s  nonml(;ratlng  bubblp.  Actually,  however,  this  model  test  as 
well  as  npproxiniatlcn  (d}  are  much  too  crude  to  allow  any  quantitative  conclusions. 

Ustni’  the  concept  of  the  average  bubble  rise  mentioned  below.  It  can  be  estimated  that  the 
third  bubblp  iiiinimum  occurred  at  a  dep*h  of  about  TOO  ft.  IX  the  bubble  would  continue  to 
onclllatc  (which  It  does  tf  there  Is  still  uncondensed  vapor),  it  would  have  to  run  through  sev¬ 
eral  mure  cycles  before  it  would  reach  the  water  surface.  At  each  minimum  more  vapor  would 
be  condensed,  and  It  seems  unlikely  that  any  substantial  amount  of  vapor  would  be  left  when  the 
bubble  reached  the  water  surface. 

This  may  seem  lo  be  in  contradiction  to  tne  observed  surface  phenomena  which  began  about 
to  Her  after  the  explosion  and  which  nhowed  quite  a  resemblance  to  the  “breakthrough"  of  a  gas 
bubble.”  These  phenomena,  however,  are  not  necessarily  caused  by  the  venting  of  a  steam 
bubble.  At  the  moments  of  the  first  and  second  bubble  maxima,  the  mass  of  water  near  the 
bubble  acquires  an  upward  momentum  which  produces  a  flow  directed  upward.  This  momentum 
will  be  conserved,  and  the  flow  win  continue  even  when  all  the  steam  is  condensed.  The  violent 
upwclllni'  of  thin  movini;  mass  of  water  probably  produced  the  plumes  and  the  surface  phenom¬ 
ena  of  Wigwam. 

Figure  3.3  shows  that  all  bubble  maxima  and  minima  lie  approximately  on  a  straight  line. 
This  means  that  the  bubble  rises  with  the  same  average  velocity  In  each  cycle.  Condensation 
would  reduce  the  periods  of  the  osclllatton,  but  the  bubble  maxima  and  minima  would  be  ap- 
proxinmtcly  on  the  same  line,  and  the  bubble  or  the  water  surrounding  it  would  retch  the 
water  surface  about  the  same  ttmc  as  the  gas  bubble.  Hence  the  fact  that  the  plume  formation 
began  at  the  time  predicted  (or  a  TNT  explosion  does  not  necessarily  indicate  a  "TNT-llke" 
bubble  behavior. 

If  all  the  vapor  were  condensed,  the  water  would  move  upward  at  a  constant  rate,  namely, 
the  aboic-mcntloncd  average  velocity,  which  Is  found  to  be  216  ft/aee.  Disregarding  air  drag, 
the  water  would  rise  72S  ft  above  the  water  surface,  which  fa  much  less  than  the  observed 
14n0-f(  maximum  plume  height.”  The  latter  corresponds  to  a  velocity  of  SOO  (t/sec,  which  Is 
40  per  cent  higher  than  'he  average  rate  of  rise.  This  sgsin  ts  not  a  proof  that  the  plumes  are 
driven  up  by  expanding  gases  or  vapor,  because  some  portions  of  the  moving  mass  of  water 
ean  have  higher  veloeiiles  than  its  rate  of  rise.  For  Instance,  a  suitable  hydrodynamic  model 
of  such  a  moving  mass  of  water  Is  a  vortex  ring.  (Vortex  rings  have  been  actually  observed 
with  high  explosives  fired  at  great  depth.)  The  total  klr-ttlv  energy  of  eucli  a  vortex  ring  ie 
larger  than  Its  translational  kinetic  energy.  When  reaching  the  water  eurface,  some  of  ite 
particles  will  rise  much  higher  than  the  average  calculated  from  the  translations!  energy. 

These  considerations  show  that  the  phenomena  observed  In  Wigwam  are  no*,  tn  contradic¬ 
tion  to  the  pt^ssibllliv  that  all  the  steam  Is  conoensco  in  the  later  bubble  oacillations.  Complete 
condensation,  however,  must  not  be  expected.  For  Instance,  all  the  gasea  dissolved  In  the 
evaporated  water  will  remain,  and  combined  with  these  will  be  a  certain  amount  of  water  vapor. 
This  Is  clear  Iv  visible  In  the  model  tests,  but  an  estimate  of  Its  .magnitude  Is  difficult  and  has 
not  been  attempted.  For  practical  purposes  this  smount  is  probably  negligibly  small,  and  it 
seems  safe  to  state  that  cs«cnttally  all  the  vapor  was  condensed  in  Wigwam  before  the  bubble 
reached  the  surface. 


3.6  SUMMARY 

The  babble  period  of  an  atomic  explosion  of  30-kt  yield  at  h  depth  of  2000  f*.  under  water 
was  calculated  to  be  2.88  sec.  A  24.3S-kl  TNT  charge  detonstel  at  Ihe  same  depth  would  have 
the  same  period. 

The  measured  bubble  period  in  Wigwam  Is  2.87  sec.  In  view  of  the  spproximatlons  made 
In  the  calculations,  such  good  agreement  was  not  expected. 

At  the  bubble  maximum,  19.9  kt  of  water  are  In  the  vqior  form.  The  liquid  producing  this 
would  make  a  sphere  of  the  same  site  as  a  30-kt  TNT  iqiheri. 

Data  lor  a  nngratini’  bubble  of  a  24,35-kt  TNT  explosion  arc  presented  Ip  Sec.  3,5  and  are 
lilu.stratcd  in  Fig.  3.3.  The  i;ax  bubble  reaches  the  water  surface  about  11  sec  after  the  detona- 
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tiori.  Comparison  with  model  tests  using  electric  sparks  as  energy  sources  Indicates  that  the 
steam  bubbles  produced  by  atomic  explosions  must  behave  somewhat  differently.  There  will  be 
condensation  of  vapor  at  the  bubble  minima.  Very  little  vapor  Is  expected  to  reach  the  water 
surface.  The  plumes  observed  In  Wigwam  about  iO  sec  alter  the  ejq^tlosion  are  probably  pro¬ 
duced  bv  the  violent  upweliing  of  water  which  originally  surrounded  the  bubble.  This  water 
keeps  on  moving  upward  even  when  all  steam  is  condenaed. 
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APPEVDIX  A 


GLOSSARY  OF  SYMBOLS 


a,, >2  Abbreviations  defined  by  Eqs.  1.38  and  1.36 

Ai  'Uh  ciiefficioit  in  the  expanaiun  uf  llie  reduced  pressure  aboul  the  center 
A;^t  Maximum  hnhhlp  radius 

B|  Uh  coefficient  in  the  expansion  of  the  reduced  density  about  the  center 
B  Depth  of  bottom  below  the  center  of  the  ertpiouion 

C|  Uh  coefficient  in  the  expansion  of  the  reduced  velocity  about  the  center 
C  Integration  constant  in  Eq.  2.4:a 

r  Velocity  of  sound;  parameter  in  TFD  theory;  constant  in  Eq.  2.32 

<\  Velocity  of  sound  1”  the  undisturbed  medium 

<•"  ^  Ideal-gas  heat  capacity  at  constant  volume  tor  the  Uh  constituent 

( °  Ideal-gas  heat  capacity  at  constant  orcssurc  for  the  Uh  constituent 

Cji/cj;  Average  heat  rapacities  defined  by  Eq.  1.16 
c  2.71828 . . . ;  elcctron'r  charge 

D  Depth  of  explosion  below  water  surface 

E  Internal  energy  per  unit  mass 

E|  Iriternai  energy  per  unit  mass  directly  behind  the  shock  front 

E«  Interna!  energy  oer  unit  nioss  in  the  undisturbed  medium 

Esw  Sho  -k-vt-avc  energy  flux 

f  Sihape  factor 

gl  Statislica'  weight  uf  the  Uh  excited  state  of  an  atom,  ion,  or  molecule 

G|  Decay  facttir  uf  the  .sh'-ck-wave  peak  density,  Eq.  2.11 

h  Dissipated  enthalpy  increment;  Planck’s  constant 

Hj  Decay  factor  of  the  shoex-wave  peak  pressure,  Eq.  2.1& 

J  Reduced  internal  energy 

Ji  Rcciuctd  internal  energy  directly  behind  the  shock  front 

j’  Reduced  internal  energy  of  an  ideal  gas 

k|  Covolume  factor  for  the  Uh  constituent 

kj  CuvoUimc  factor  for  a  component  when  in  its  ground  state 

k  Total  covolume  factor 

k'  Boltzmann's  constant 

K|  Equiliorium  ronstanl  for  the  jf/i  chemlcai  reaction 

K|,_  Ideal-gas  equilibrium  constant  for  the  ]f/t  chemical  reaction 

in  Natural  logarithm 

log  logarithm  to  the  base  10 

l.i  Decay  factor  for  the  shock -wave  peak  particle  velocity,  Eq.  2.12 

r.i  --  3  fr’  -  J);  particle  mass;  total  mass  ol  the  sphere  encompassed  by  the  shock  front 
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M(  Molecular  wcljjlit  oJ  water 

M  Number  of  hydrogen  nuclei  per  cell  In  the  TFD  theory 

ni  Number  of  moles  of  the  Uh  constituent  in  grrtns  of  the  mixtore 

nt>  Total  number  of  moles  In  M«  grains 

Nav  Avogadro’s  number 

p  Excess  pA'essure 

Pj  Shock-wave  peak  pressure 

P^  Thermal  pressure 

Pj  Absolute  pressure  In  the  undisturbed  medium 
P(j  Coe'flclents  in  Eqs.  2.46  and  ?.47 

Pm  Pressure  In  the  steam  bubble  at  the  moment  of  the  maximum  radius 

Ps  Pressure  of  the  saturated  liquid 

Q  Energy  yield  In  kllotons  of  TNT;  partition  function 

Qdlaa  Dissipated  energy 

Wd  nubble  energy 

r|]  Coefficients  In  Eq.  1.24 

r  Radius 

T]  Radius  of  shock  front 

r|  Radius  of  the  layer  of  particles  which  forms  the  bubble  Interface  at  Its  maximum 
expansion 

rg  Reference  radius 

R  Qas  constant 

(R  Cell  radius  In  TFD  theory 

S  Entropy  per  unit  maes 

&t  Entropy  per  unit  mass  directly  behind  the  shock  front 
Sg  Entropy  per  un't  mass  in  the  undisturbed  medium 

Sam  Entropy  per  unit  mass  at  the  bubble  maximum 

Scr  Entropy  per  unit  mass  at  the  critical  point 

t  Time 

T  Temperature;  first  bubble  -erlod 

Tg  Free-water  bubble  perloe 

u  Particle  velocity 

U|  Particle  vetociiy  directly  behind  the  shock  front 

ug  Particle  velocity  In  the  undisturbed  medium 
U  Propagation  velocity  of  the  shock  front 
V  Specific  volume 

Vg  Specific  volume  of  the  saturated  liquid 

Vg  Ei'fecttve  atomic  or  ionic  volume  of  a  component  In  its  ground  state 

V|  Effective  atomic  or  Ionic  volume  of  a  component  when  In  Its  Ifh excited  state 

W|  Coefficients  In  £q.  2.4P 

W  Energy  yield  In  pounds  of  TNT 

X  Imperfection  factor  In  the  HKW  equation  of  state,  Xq.  1.17;  r.-.'ucsd  radius  In  TFD  and 

Snay-Matthixs  theories 

*  =  (P|  -  Pg)/pg 

Vi  Number  of  moles  of  hydre^en  stoms  In  Mg  grams  of  the  mixture 

y,  Numoer  of  motes  of  oxygen  atoms  In  Mg  grams  of  the  mixture 

yg  Number  of  moles  of  free  electrons  In  Mg  grams  of  the  mlxturr. 

Z  Charge  of  the  oxygen  nucleus 

Zi  Valency  ol  the  Uh  component 

Zj  Valency  of  the  Uh  reaction 

Zg  Total  hydrostatic  head  In  feet 

a  T.me  factor;  parameter  In  HK.W  equation  of  state 

ayi  Factor  which  accounts  for  the  Internal  energy  within  the  bubble 

Logarithmic  decay  fac'or  lor  the  reduced  total  energy;  parameter  In  uicW  equation  ol 
rtatc 
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(?o  Exchange  correction 

0s  Logarlthinlc  decay  factor  lor  the  reduced  shock-wave  energy 

V  Isentroplc  exponent,  Eq.  1.9 

r,  Isentroplc  exponent  directly  l  ehlnd  the  shock  front 
Isentroplc  exponent  for  an  Ideal  eas 

mp  Logarlti.mlc  derivative  of  the  peak  pressure  with  respect  to  the  peak  density 

€|  Excitation  energy  of  the  1/ft  excited  state  of  an  atom,  Ion,  or  molecule 

i)t  Reduced  total  energy,  Eq.  2.3 

t)S  Reduced  shock-wave  energy 

{  Reduced  radius 

VI  Stoichiometric  coefficient  of  the  1/ft  component  In  a  chemical  reaction 

s  3.141585... 

,0  Density 

Pi  Density  directly  behind  the  shock  front 
Pg  Deiiaiiy  of  itie  uiiciBiui'tieu  uiuiiluui 

p  lliomas-Ferml  unit 

a  Coel'iclent  In  Eq.  1.34 

*  Electrostatic  correction,  Eq.  1,23 

(p  Reduced  particle  velocity 

<Pi  nth  derivative  of  <p  with  respect  to  evaluated  ai  (  =  1 
X  Reduced  density 

Ip  Reduced  pressure 

♦  TFD  potential 

0  Time  constant 
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The  Surgeon  •jeneral,  D/A,  Washington  25,  D.  C.  ATTN:  Chlel,  X&D  Olviilon  7 

Chief  Chemical  Officer,  D/A,  Waihlngton  25,  D.  C.  B-3 

The  Quarterrnatter  General,  D/A,  Washington  25,  D.  C.  ATTN:  Research  and  Development  Dlv.  10 

Chief  of  Engineers,  D/A,  Washington  26,  D,  C.  ATTN:  ENCN*  11-16 

Chief  of  Transporutlon,  Mlltury  Planning  and  Intelligence  Div»  Washington  25,  D.  C.  16 

Commanding  General,  ContlneiiUl  Army  Command,  Ft.  Monroe,  Va.  17-10 

Presldei.t,  Board  51,  lleadduarters,  Contlnenul  Army  Command,  Ft.  Sill,  Okla.  20 

President,  Board  *2,  Headipianeta,  Continental  Army  Commend,  Ft.  Knox,  Ky.  21 

President,  Board  #4,  Headquarters,  Continental  Army  Command,  Ft.  Bliss,  Tex.  27 

Commanding  General,  U,  S.  Army  Caribbean,  Ft.  Amador,  C.  Z.  ATTNt  Cml.  Off.  23 

Commander-m-Chlef,  Fat  East  Coritmand,  APOiOO,  Sen  Francisco,  attv-  s/-.,rc  t  ■>  24-15 

Commanding  General.  U.  S.  Army  Europe,  APO  403,  New  York,  N.  Y.  ATTNt  OfOT  Dlv„ 

Combat  Dev.  Bt.  26-27 

Commanding  General,  U.  S.  Army  Pacific,  APO  056,  .San  Francisco,  Calif.  ATTNt  Cml.  Off.  26-2B 

Commandant,  Command  and  Gecetal  Staff  College,  Ft.  Leavenwonh,  Kans.  ATTNi  ALLLSCAS)  30-31 

Commandant.  The  Artillery  and  Guided  Missile  School,  Ft.  SUl,  OHa.  32 

Secretary,  The  Antiaircraft  Artillery  and  Gulden  Missile  School,  Ft.  Blla,  Tex.  ATTN;  MaJ  Gregg  D. 

Breitegan,  Dept,  of  Tactics  and  Combined  Aims  33 

Commanding  General,  Army  Medical  Service  School,  Brooke  Army  Medical  Center,  Ft.  Sam 
Houston,  Tex.  34 

Director,  Special  Weapons  Development  Olflce.  Headquarters,  CONAIC,  Ft.  Bliss,  Tex.  ATTN: 

Caps  T,  E.  Skinner  36 

Commsndant.  Walter  Reed  Army  Institute  of  Research,  Walter  Reed  Anny  Medical  Center, 

Washington  25,  D.  C.  36 

Superintendent,  11.  S.  MlUtaty  Academy,  West  Point.  N.  Y,  ATTN:  Prof,  of  Ordnance  37 

Commandant,  Chemical  Corps  Schcol.  Chemical  Corps  Training  Command,  Ft.  McClellan,  Ala.  38 

Commanding  General,  Research  and  Engineering  Command,  Army  Chemical  Center,  Md.  ATTN: 

Deputy  for  RW  and  Non-Toxic  Moerlat  39 

Commanding  General,  The  Et.glneer  (.enter.  Ft.  selvou,  Va.  ATTN:  Aut.  Commanda.tt,  Engince, 

School  10-42 

Commanding  Oftlcer,  Engineer  Research  and  Development  Laboratory,  Ft.  Belvoir,  Va.  ATTN; 

Chief,  Trch.ni'-al  Intelligence  Branch  43 
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(:■  Ofttcrr,  rit»ilrny  Ar:rn*l,  Dover,  N.  J,  ATTNi  OKDH'TK 

CojKiiundiiig  Officer,  Army  Modlcut  Reiclfcli  I.tborttory.  Pi.  Knox,  Ky. 

Coiirmindliig  Officer,  ClicmlcAl  Ccrpi  Chemli'.Al  end  lUdlological  Uborttoiy,  Army  Chomlcol 
Center,  Md.  ATTNi  Tech.  Ubr«'y, 

Commending  Officer,  Tramportetlon  RAD  Stillon,  Pu  EutUi,  V<. 

Director,  Teclinical  Documenu  Center,  Ev«m  llgnil  UboMUty.  lolmor,  N. ). 

Director,  Wtterwayi  Experiment  Sution,  PO  Inx  631,  Vlekibut|,  MfM.  ATTNi  Ubtaiy 
Director,  Armed  Forcei  liuiitutc  of  Ptt'tcl''gy,  Welter  heed  Army  Mndlcel  Conior,  6638  160i  Street, 
N.W.,  Wiihington  26,  D.  C. 

Director,  Opeietloni  Reicetch  Office,  lohhi  HopUiu  Urtlvcnity,  3100  Connecticut  Avc«  Chevy  Chaie, 
Md.,  Wtihlngton  IS,  D.  C. 

CommAndlng  Oenerel,  QuertermiKet  Rcecetcit  end  Development  Commend,  OuiitenuMer  RMeercb 
end  Development  Center,  Natick,  Mata,  ATTNi  CBR  Ualaon  Officer 

NAVTf  ACTIVITIES 

Chief  of  Naval  (.tperctloni,  D/N,  Waihlngton  36,  D.  C.  ATTNi  OP>38 
Chief  of  Naval  Opcratloiti,  0/N,  Wuhlngton  38,  D.  C.  ATTNi  OP-OSEC 
Director  of  Naval  Intelligence,  D/N,  Wathlngto:.  26,  D.  C.  ATTNi  OP>033V 
Chief,  Dureau'of  Medicine  tnd  Surgery,  D/N,  Waihlngton  36,  D.  C.  ATTNi  Speclel  Wetpon 
Dcfcnie  Div. 

Chief.  Bureau  of  Ordnance,  0/N,  Waihlngton  25,  D.  C. 

Chief.  Bureau  ofShipi,  0/N,  Waihlngton  36,  D.  C.  ATTNi  Code  348 
Chief,  Bureau  of  Yardi  and  Docki,  D/N,  Waihlngton  26,  D.  C.  ATTNi  D<440 
Cldef,  Bureau  of  Suppllei  and  Accounti,  D/N,  Waihlngton  25,  D.  C. 

Chief,  Bureau  of  Aeronautlci,  D/N,  Waihlngton  25,  0.  C. 

Chief  of  Naval  Reicarch,  Department  of  the  Navy,  Waihlngton  25,  D.  C.  ATTNi  Code  811 
Commander-In-Chief,  U.  S.  Pacific  Fleet,  Fleet  Poit  Office,  San  Franclico,  CtUf. 
Commander-In-Chief,  U.  S.  Atlantic  Fleet,  U,  S.  Navel  Bale,  Norfolk  11,  Va, 

Commandant,  U.  S.  Marina  Corpi,  Waihlngton  36,  D,  C.  ATTNi  Code  AOSH 
riciident,  U.  S.  Naval  War  College,  Newport,  R.  L 
Superintendent,  U.  S.  Naval  Poitgiaduate  School,  Monterey,  Calif. 

Commanding  Officer,  U.  S.  Naval  Schooli  Command,  U.  5.  Naval  Sutiou.  Tteaaure  Uland, 

San  Pranclico,  Calif, 

Commanding  Officer,  U.  S.  Fleet  Training  Center,  Naval  Bate,  Norfolk  11,  Va.  ATTNi  Special 
Wcaponi  School 

Commanding  Officer,  U.  S.  Fleet  Training  Center,  Naval  Station,  Ban  Diego  36,  CalU,  ATTNi 
(SPWP  School) 

Commanding  Officer,  1.1,  S.  Naval  Damage  Control  Training  Center,  Naval  Bate,  PMlidelphla  13. 

Pa.  ATTN:  ABC  Defenie  Coune 

Commanding  Officer,  U.  S.  Naval  Unit,  Chemical  Corpi  School.  Army  Otcmlcal  Training  Center, 

Ft.  kIcCiellan,  Ala. 

Commander,  U.  S.  Naval  Ordnance  Laboratory,  Silver  Spring  13,  Md.  ATTNi  EE 

Commander,  U,  S.  Naval  Oiunaucc  Laboratory,  Sliver  Spring  19,  Md.  ATTN;  EH 

Commander,  U.  S.  Naval  Ordnance  Laboratory,  Silver  Spring  19,  Md,  ATTNt  R 

Commander,  U.  S.  Naval  Ordnance  Ten  Station,  Inyokern,  China  Lake,  Calif. 

0(flccr-in -Charge,  U.  S.  Naval  ClvU  Engineering  Rea.  and  Evaluation  Labw  l>.  S.  Naval  CaiMiuctloa 
Banalion  Center,  Port  lluencme,  Calif.  ATTNi  Code  363 
Commanding  Officer,  U,  S.  Naval  Medical  RcHarch  Inn.,  National  Naval  Medical  Center,  Betheida 
14,  Md. 

Director,  Naval  Air  Expertmenul  Sution,  Alt  Material  Cantar,  U.  I.  Naval  Bam,  PUladtlphla,  Fa. 
Director,  U,  S.  Naval  Reiearch  Laboratory,  Waihlngton  38,  D.  C.  ATTNt  Mn.  KatheilBe  H.  Can 
Commanding  Cificer  and  Dlrecioi,  U.  S.  Navy  Elcctronict  Laboiatoty,  San  Diago  58,  Calif. 

ATTNt  Code  4223 

Commanding  Officer,  U.  S.  Naval  Radiological  Dcfanar  Laboratory,  tan  ftanelMo  B4.  Calif. 

ATTNi  Technical  Information  Dlvlilon 

Commanding  Officer  and  Director.  David  W.  Taylor  Modal  Batin,  WaMngton  3,  D.  C.  ATTNi 
Library 

Commander,  U,  8.  Naval  Air  Development  Center,  lohmvllle,  Pa. 
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C'  'T.nii.n’Jlim;  Ofitcrr,  ritailrny  Af:.fn»l,  Dover,  N.  I.  ATTNi  OIU>M«TK 
(.'•'iKiiunditig  Officer,  Army  Mcdlcul  ReKtrcIi  l.tborctory.  Pi.  Knox,  Ky. 

Commtndliig  Officer,  ChcmlcAl  Cerpi  Chemli'.xl  tnd  lUdiologicxl  Lobotttoiy,  Army  Chomicxl 
Center,  Md.  ATTNi  'I'ccli.  Ubra-y, 

Commanding  Officer,  Tranrportatlr^rn  HAD  Station,  Pu  Euttii,  Va. 

Director,  Teclinicat  Documenta  Center,  Evana  Signal  Laboratory,  Mmar,  N.  i. 

Director,  Waterwayi  Experiment  Sution,  PO  Inx  631,  Vlehtbutg,  MiH.  ATTNt  Ubraiy 
Director,  Armed  Parcel  liuiitute  of  Pat'tcl''gy,  Walter  keod  Army  Mndlcal  Conur,  6636  16rJi  Strtet, 
N.W.,  Waihlngton  26,  D.  C. 

Director,  Opeiationa  Reicatch  Office,  lohni  Hopklna  Unlvcnity,  3100  ConnoeUcut  Avo«  Cbovy  duac, 
Md.,  Waihington  15,  D.  C. 

CommAndlng  Ceneral,  Quartermaner  Reacarcli  and  Oevelopmont  Command,  Quaitatmaator  Raaearcb 
and  Development  Center,  Natick,  Maaa.  ATTNt  CBR  Lialaon  Offlcor 

NAV'Y  ACTIVITIES 

Chief  of  Naval  Operetlona,  D/N,  Waihlngton  36,  D.  C.  ATTNt  OP>38 
Chief  of  Naval  Operation!,  0/N,  Waahlngtoil  36,  D.  C.  ATTNi  OP-03EC 
Director  of  Naval  Intelligence,  D/N,  Waihlngto:.  36,  D.  C.  ATTNi  OP-039V 
Chief,  Dureau'of  Medicine  and  Surgery,  D/N,  Waihlngton  36,  D.  C.  ATTNt  Specif)  Weaponi 
Defenie  Dlv, 

Chief,  Bureau  of  Ordnance,  D/N,  Waihlngton  26,  D.  C. 

Chief,  Bureau  of  Shlpi,  D/N,  Waihlngton  36,  D.  C.  ATTNt  Code  346 
Chief,  Bureau  of  Yardi  and  Docki,  D/N,  Waihlngton  20,  D.  C.  ATTNt  D-440 
Cliief,  Bureau  of  Suppliei  and  Accounti,  D/N,  Waihlngton  26,  D.  C, 

Chief,  Bureau  of  Aeronautlui,  D/N,  Waihlngton  26,  D.  C. 

Chief  of  Naval  Ikicarch,  Department  of  the  Navy,  Waihlngton  26,  D.  C.  ATTNt  Code  811 
Commander-In-Chief,  U.  S.  Pacific  Reet,  Fleet  Poit  Office,  San  Franctaco,  Calif, 
Commandet-in-Chief,  U.  S.  Atlantic  Fleet,  U,  S.  Naval  Baae,  Norfolk  11.  Va. 

Commandant,  If.  S,  Marine  Corpi,  Waihlngton  36,  D.  C.  ATTNt  Code  A03H 
rreiident.  U.  S,  Naval  War  College,  Newport,  R.  L 
Superintendent,  U.  S.  Naval  Pottgraduate  School,  Monterey,  Calif. 

Commanding  Officer,  U.  S.  Naval  Schoola  Command,  U,  S.  Naval  Sutiou,  Tieaaur*  laland, 

San  Ptanclico,  Calif, 

Commanding  Officer,  U.  S.  Fleet  Training  Center,  Naval  Bate,  Norfolk  11,  Va.  ATTNt  Special 
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Coitirnandlng  Offlcei,  U.  S.  Fleet  Training  Center,  Naval  Station,  Ban  Diego  36,  CalU,  ATTNt 
(SPWP  School) 

Commanding  Officer,  II.  S.  Naval  Damage  Conrial  Training  Center,  Naval  Bale,  FUlidelphla  13. 

Pa.  ATTN:  ABC  Defenie  Coune 

Commanding  Officer,  U,  S.  Naval  Unit,  Cticmlcal  Corpi  Scnool.  Army  Gtcmical  TtaUlng  Center, 

Ft.  ktcCiellan,  Ala. 

Commander,  U.  S.  Naval  Ordnance  Laboratory,  Silver  Spring  19,  Md.  ATTNt  BE 

Commander,  U,  S.  Naval  Oiduauce  Laboratory,  Silver  Spilng  19,  Md.  ATTN:  EH 

Commander,  U.  S.  Naval  Ordnance  Laboratory,  Silver  Spring  IB,  Md.  ATTNt  R 

Commander,  U.  S.  Naval  Ordnance  Ten  Station,  inyoketn,  China  Lake,  Calif. 

Officcr-ln -Charge,  U.  S.  Naval  Civil  Engineering  Rei.  and  Evaluation  Labn  D.  S.  Naval  Coumuctioa 
Battalion  Center,  Port  lluencmc,  Calif.  ATTNt  Code  363 
Commanding  Officer,  U,  S.  Naval  Medical  RcMarch  Inii.,  National  Naval  Medical  Cemer,  Bciheada 
14,  Md. 

Director,  Naval  Air  Expertmenut  Station,  Alt  Material  Center,  U.  3.  Naval  Baae.  PUladelphln,  Fa. 
Director,  U,  S.  Naval  Reicarch  Laboratory,  Waihlngnn  36,  D.  C.  ATTNt  Mn.  KatbeiiM  H.  Caaa 
Commanding  Cificer  and  Director,  U.  S.  Navy  EIcctronict  Laboratory,  San  Diego  83,  Calif. 
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Commanding  Officer,  U.  S.  Naval  Radiological  Dcfcna*  Laboratory,  San  ftanelKO  34.  Calif. 

ATTNi  Technical  Information  Dlvlilon 

Commanding  Officer  and  Director.  David  W.  Taylor  Model  Bciin,  Wadilngten  3,  D.  C.  ATTNi 
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Commander,  U,  S.  Naval  Air  Development  Center,  lohnavllle,  Fa. 
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Commander,  Operational  Developmeut  Botce,  United  Sutea  Atlantic  rUet,  Unllad  Siatoc  Kivai 
Balk,  Norfolk  11,  Va.  BBS 
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